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Resumen: 

Con el propósito de suministrar un marco de referencia adecuado para el desarrollo de los modelos 
matemáticos que son el objetivo principal de esta investigación, en la primera parte de este artículo, 
se resumen los trabajos de campo llevados a cabo en las aguas costeras del Río de la Plata en la 
región de Montevideo y en la Bahía de Montevideo, entre 1991 y 2001. Se comienza por los resul- 
tados de un estudio de la dinámica y renovación del agua en la bahía de Montevideo que involucró 
el empleo de trazadores radioactivos y fluorescentes para marcar el agua en el interior de la bahía. 
Durante ese estudio, la detección de una señal radioactividad artificial debida al tritio a 15 km al 
este de la boca de la bahía, condujo a conjeturar la existencia de un proceso de difusión de rompien- 
tes en la franja de aguas costeras adyacentes, situadas al este de dicha boca. Para distinguir entre los 
aportes del emisor sub-acuático de Punta Brava y los aportes de la bahía se realizaron estudios de la 
influencia de la bahía sobre las aguas costeras en simultáneo con un estudio de la dinámica de la 
pluma de descarga del emisor sub-acuático. Durante estos trabajos se descubrió un efecto de trampa 
significativo en la playa Pocitos. Se investigó este efecto empleando trazadores fluorescentes para 
simular un evento localizado de contaminación. A continuación se efectuó un estudio de dinámica 
de sedimentos en la bahía de Montevideo y en las aguas costeras adyacentes, empleando trazadores 
radioactivos. Además se obtuvo información relacionada con procesos de dragado empleando son- 
das de densidad. Se determinó la capacidad de autodepuración de las aguas costeras en la región de 
Punta Lobos, donde se ha previsto ubicar el segundo emisor sub-acuático de la ciudad de Montevi- 
deo, empleando un método de doble trazador. 

En la segunda parte del artículo se revisan y discuten los modelos matemáticos (modelos analíticos) 
construidos y utilizados para diseñar y para interpretar los resultados de los trabajos de campo 
resumidos en la primera parte. Se comienza por los modelos de compartimentos construidos para 
una caracterización paramétrica de la cinética de renovación del agua en la bahía de Montevideo, 
incluyendo un modelo clásico de dos compartimentos y un modelo reciente de tres compartimentos 
basado en una cinética fractal. Se resumen los modelos matemáticos utilizados en los trabajos de 
campo en aguas costeras, ya sea con fines de pronóstico para protección radiológica y diseño de 
experimentos, ya sea con fines de diagnóstico para la caracterización paramétrica de los procesos de 
transporte a partir de los resultados experimentales. Se revisa el modelo matemático construido para 
la determinación del Too en las aguas costeras vecinas a Punta Lobos. 


Palabras clave: modelos matemáticos, trazadores, efecto de escala, cinética fractal, difusión anó- 
mala, capacidad de autodepuración. 


Abstract: 
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With the purpose of furnishing a suitable reference framework to develop the mathematical models 
that are the main objective of this research, the first part of the article summarizes the fieldwork 
carried out in the coastal waters of the Río de la Plata in the Montevideo region and in the Bay of 
Montevideo between 1991 and 2001. We begin with the results of a study of the dynamics and 
renewal of water in the bay of Montevideo that involved the use of radioactive and fluorescent 
tracers injected in the water inside the bay. During this study, the detection of artificial radioactivity 
due to tritium, 15 km east from the bay?s mouth, led to the conclusion that there is a process of 
breakers diffusion in the strip of adjacent coastal waters located to the east of said mouth. In order 
to distinguish between the contributions of the Punta Brava sub-aquatic sewage outfall and the 
contributions of the contamination of the bay, studies of the influence of the bay on the coastal 
waters were carried out simultaneously with a study of the dynamics of the sewage outfall discharge 
plume. During these works a significant trap effect was discovered at Pocitos beach. This effect was 
investigated using fluorescent tracers to simulate a localized contamination event. A study of sedi- 
ment dynamics was carried out in the Bay of Montevideo and adjacent coastal waters using radioac- 
tive tracers. In addition, information related to dredging processes was obtained using density 
probes. The self-purification capacity of the coastal waters in the Punta Lobos region (where it is 
planned to locate the second sub-aquatic emitter of the city of Montevideo) was determined using a 
double tracer method. 

In the second part of the article we review and discuss the mathematical models (analytical models) 
constructed and used to design and interpret the results of the field work summarized in the first 
part. We start with the compartment models constructed for a parametric characterization of the 
water renewal kinetics in the Bay of Montevideo, including a classical two compartments model 
and s recent model three compartments model based on fractal kinetics. We summarize the mathe- 
matical models used field work in coastal waters, either for forecasting purposes (for radiological 
protection and design of experiments) or for diagnostic purposes (for parametric characterization of 
transport processes from experimental results). The mathematical model constructed for the deter- 
mination of Too in the coastal waters adjacent to Punta Lobos is reviewed. 


Key words: mathematical models, tracers, scale effect, fractal kinetics, anomalous diffusion, self- 
purification capacity. 


Introducción! 


Las aguas costeras cuyo estudio motiva el presente trabajo se encuentran ubicadas en la 
margen izquierda del Río de la Plata, en el inicio de su región externa. La marea, las des- 
cargas fluviales y los vientos presentan efectos variables sobre la dinámica del agua en el 
Río de la Plata en general y de las aguas costeras de Montevideo en particular, pero de un 


| Una parte del contenido de este artículo se basa en una conferencia plenaria presentada por el 
autor en el octavo encuentro del International Center for Earth Sciences, realizado en Mar del Plata, 
Argentina, en 2012 [1]. Exceptuando la construcción y aplicación de modelos de mecánica de flui- 
dos computacional (CFD) en la bahía y en las aguas costeras próximas, la totalidad de los trabajos 
de campo, de laboratorio y de modelado matemático que se relatan fueron coordinados por el autor 
y llevados a cabo con la participación de diversas instituciones públicas y privadas de Uruguay 
(Ministerios, Universidades, Fuerzas Armadas, Administración Nacional de Puertos, Intendencia 
Municipal de Montevideo, ADES). Además, en la mayor parte de los trabajos participaron expertos 
del Organismo Internacional de energía Atómica de NNUU a través de proyectos nacionales de 
cooperación que se mencionan en el desarrollo del artículo. Los modelos originales que aparecen en 
el presente trabajo fueron construidos por el autor. 
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mismo orden de magnitud. Como consecuencia, tanto desde un punto de vista hidrológico 
como ambiental, las aguas costeras de Montevideo configuran un sistema complejo y 
cambiante que complica el control de la calidad del agua de las playas. 


El relieve de Montevideo presenta ondulaciones, con algunas alturas suaves. Por ello está 
cruzado por numerosos y por lo general pequeños cursos de agua. Esta red hidrográfica 
tiene mucha importancia porque ha condicionado fuertemente el desarrollo del sanea- 
miento y la disposición final de las aguas servidas de la ciudad. Con excepción de una 
franja costera de un ancho de unos 3 km, que vierte directamente hacia el Río de La Plata, 
el resto de la ciudad drena hacia tres arroyos. El Pantanoso, con 16 km de extensión total 
y una cuenca de 66.4 km? y el Miguelete, de 21,5 km de extensión y una cuenca de 113 
km?, desaguan en la Bahía de Montevideo. El arroyo Carrasco, de 14 km de extensión y 
una cuenca de 213 km? recoge aguas de extensas áreas pobladas del este y noreste de 
Montevideo y las vuelca en el Río de La Plata. 

Estos tres cursos de agua están contaminados y degradados por numerosas descargas de 
aguas servidas domésticas y residuales industriales, así como también por basurales. La 
existencia en la zona de decenas de industrias, cuyos efluentes llegaban hasta hace pocos 
años, directa o indirectamente al curso de agua constituye uno de los factores determinan- 
tes de la contaminación de la Bahía de Montevideo (Figura 1). 


zo 


Ss 


VIN3LSI 


AWNLIV NOIDWNLIS 
30 S 


OLNIINVINVS 


OJUIAJLNON 301 OLNINVLIVIJO 


ZA men severa 5, RIO DE LA PLATA 


Figura 1. Esquema del saneamiento de Montevideo en el período 1991-1995. 
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En la figura se indican los arroyos Carrasco (desembocando en el Río de la Plata al Este), 
Pantanoso y Miguelete desembocando en la bahía. En Punta Carretas se encuentra el 
actual emisor sub-acuático, mientras que en Punta Lobos, en el extremo Oeste de la boca 
de la bahía, se ubicará un segundo emisor sub-acuático de la ciudad. 


En 1990 comienza a funcionar el emisor subacuático, que recolecta las aguas y las vuelca 
río adentro, a una distancia de 2.200 metros al sur de Punta Brava (Punta Carretas). 

La calidad de las aguas de las playas de Montevideo mejora con la entrada en servicio del 
emisor subacuático, pero menos de lo esperado. Una explicación podía ser que el emisario 
no funcionara de acuerdo a lo previsto. Pero si lo hacía bien, entonces debía haber otra 
razón para la contaminación residual. 

Una posibilidad para explicar esa contaminación remanente era que se debiera a ciertos 
vertidos directos de aguas residuales en las cercanías de la orilla, en el arco de playas al 
Este de la boca de la bahía, vertidos que todavía no habían sido incorporados en un colec- 
tor perimetral que los condujera al emisor subacuático. 

Por otra parte, todavía hoy, cuando la lluvia es lo bastante copiosa, el sistema de sanea- 
miento libera las aguas residuales a través de un conjunto de aliviaderos situados en la 
costa, debido a lo cual la contaminación de las aguas costeras aumenta transitoriamente 
durante uno o dos días. 

Finalmente, otra posibilidad era un efecto contaminante debido a las aguas de la propia 
bahía de Montevideo. 

La necesidad de dilucidar el origen de la contaminación residual motivó una parte de los 
experimentos con trazadores radioactivos? considerados en el presente artículo. 


Entre 1991 y 2001 se llevaron a cabo trabajos de campo en la bahía y en las aguas coste- 
ras de la región de Montevideo, destinados a estudiar procesos de transporte de masa y 
dinámica de sedimentos. 

Para diseñar los trabajos de campo, valorar aspectos de protección radiológica de la po- 
blación y de los trabajadores ocupacionalmente expuestos, y para interpretar los resulta- 
dos obtenidos mediante una caracterización paramétrica de los procesos, se utilizaron 
modelos matemáticos analíticos. 

Algunos de estos modelos (el modelo de dos compartimentos y el modelo fractal de tres 
compartimentos para la cinética del THO en la Bahía de Montevideo, y el modelo para la 


2 Los trazadores, radioactivos o no, son una herramienta utilizada para estudiar experimentalmente 
el transporte de masa en procesos industriales, hidrológicos y ambientales. Un trazador bien 
seleccionado se integra en el medio que se desea investigar. Es estable, o si no lo es, su velocidad de 
decaimiento se conoce bien. Se puede distinguir de las partículas del medio que lo rodea por sus 
propiedades físicas o físico-químicas. Su presencia no afecta en forma significativa la dinámica 
investigada y el movimiento de las partículas de trazador resulta representativo del movimiento de 
las partículas del medio al cual se ha integrado. Entonces, estudiando cómo se distribuye 
espacialmente y cómo varía en el tiempo la concentración de trazador en el medio, se obtiene 
información sobre el comportamiento en dicho medio. 
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estimación del Too en aguas costeras estratificadas) se construyeron a la medida de las 
necesidades específicas planteadas por los experimentos, y son originales. 


PRIMERA PARTE: RESUMEN DE LOS TRABAJOS DE CAMPO 


En noviembre y diciembre de 1991 se llevó a cabo, coordinada desde la Dirección Nacio- 
nal de Tecnología Nuclear del Ministerio de Industria, Energía y Minería de Uruguay, y 
en el marco de un proyecto nacional de cooperación técnica con el OIEA, una investiga- 
ción experimental preliminar de la dinámica y renovación del agua en la Bahía de Monte- 
video [2]. La Figura 2 muestra la bahía y la región de aguas costeras adyacentes. 
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Figura 1. Situación de los principales vertidos de aguas residuales 
existentes en la zona costera de Montevideo. 


> 


Figura 2. Esquemas de descargas de aguas residuales, cloacales y pluviales, en la bahía y la costa al 
Este de la boca de la bahía en el período 1991-1995. Se indican los arroyos Miguelete y Pantanoso 
y el emisario subacuático. A la izquierda del canal de conexión del puerto con el muelle de la refi- 

nería de petróleo se observa la isla Libertad. 


Se utilizó tritio como trazador principal. Se obtuvieron resultados que sugerían que los 
contaminantes presentes en las aguas de la bahía (disueltos o integrados en sedimentos en 
suspensión) podían afectar la calidad de las aguas de las playas situadas al Este. 
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Como estos experimentos se habían diseñado para estudiar la dinámica y la renovación 
del agua en la bahía y no para estudiar un posible efecto contaminante sobre las playas, 
dada la importancia que este último proceso podría tener, se suspendió el estudio de di- 
námica de sedimentos en la bahía y se decidió estudiar la dinámica de las aguas costeras 
al Este de la boca de la bahía. 

Por ello, entre 1993 y 1996 se realizó una nueva serie de experimentos con trazadores, 
especialmente diseñados para estudiar el problema de una posible contaminación de las 
playas debida a la bahía y para obtener información sobre el funcionamiento del emisor 
subacuático de Punta Carretas. 


La primera fase de esa investigación se llevó a cabo en marzo y abril de 1993, a solicitud 
de autoridades de la Universidad de la República, utilizando tritio y fluoresceína sódica, 
también conocida como uranina. 

Simultáneamente, respondiendo a una solicitud del gobierno de Montevideo, se llevó a 
cabo un estudio experimental de la ubicación y la estabilidad de la pluma de descarga del 
emisor subacuático de Punta Brava, utilizando *%!I como trazador. Se puso a prueba la 
hipótesis según la cual dicho emisor no sería el causante de la contaminación residual 
observada en las playas de Montevideo, situadas al Este de la boca de la bahía, durante la 
temporada estival [3], [4]. 


Como puede verse en la Figura 2, inmediatamente hacia el Este de Punta Carretas (Punta 
Brava) se encuentra la Playa Pocitos. 

Los resultados de los estudios de los procesos de transporte a lo largo de la franja costera 
sugerían que esa playa presenta un considerable efecto de trampa, debido probablemente 
a la perturbación que la Punta Brava produce sobre el flujo de las aguas costeras. 

Un accidente de contaminación que se produzca en la capa superior de las aguas costeras 
adyacentes a la Playa Pocitos, si bien adquiriría un carácter marcadamente local, por eso 
mismo podría tener consecuencias graves en esa zona. 

Por este motivo en marzo, abril, mayo y junio de 1995 se investigó el efecto de trampa en 
la playa Pocitos y se determinaron los parámetros de dispersión en las aguas costeras del 
Montevideo adyacentes a Punta Brava [5]. 


Durante algunos experimentos preliminares se observaron efectos de histéresis muy mar- 
cados (la forma inicial de la mancha de trazador influía considerablemente en la posterior 
dinámica de la dispersión, hasta que la mancha de trazador superaba los 500 metros de 
diámetro). Para minimizar el efecto de histéresis se ideó un nuevo procedimiento de in- 
yección del trazador, que se empleó en los experimentos definitivos. 

Por otra parte los resultados de esos experimentos preliminares sugerían prestar atención 
a algunos detalles de la microcirculación en la capa de agua superior (particularmente en 
superficie), e intentar cuantificar un posible efecto de escala asociado al proceso de dis- 
persión en esa capa. 

Para estudiar esos detalles asociados con micro-circulación local, se efectuó una graba- 
ción de video de la evolución de una mancha de Fluoresceína Sódica, a bordo de un heli- 
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cóptero y con cámaras especiales. La digitalización de la imagen permite aplicar técnicas 
de geometría fractal a los bordes de la mancha [6]. 


En 1996 y 1997 se llevó a cabo un estudio experimental de la dinámica de sedimentos en 
la Bahía de Montevideo y en las zonas de descarga del material dragado, que había sido 
postergado luego de las campañas de 1991 para el estudio de la dinámica y renovación del 
agua. Se utilizó '"%Au como trazador para sedimentos en suspensión y !Ir para sedimen- 
tos de fondo. 

Se determinaron perfiles verticales de distribución de sedimentos, en suspensión y acu- 
mulados en la cántara de una draga de succión, empleando una sonda nuclear construida 
especialmente para este fin en el marco de un proyecto nacional de cooperación técnica 
con el OIEA [7]. 


Finalmente, nuevamente respondiendo a una solicitud del gobierno de Montevideo, en 
junio y julio de 2001 se estudió la capacidad de autodepuración de las aguas costeras de 
Montevideo en la región de Punta Lobos, en vista de la instalación de un segundo emisor 
sub-acuático. Se utilizaron Rodamina WT como trazador principal, uranina como trazador 
secundario y una población bacteriana de E. coli como contaminante microbiológico 
representativo [8]. 


Trabajos de campo 


Durante las campañas de estudio de la dinámica y renovación del agua en la Bahía de 
Montevideo se llevaron a cabo: 


-Mediciones de perfiles de temperatura y conductividad en diferentes estaciones de mues- 
treo en la bahía. 


-Medida de flujos de agua utilizando trazadores fluorescentes. 


-Marcado del compartimento noroeste de la bahía (situado por encima de una línea que 
une el espigón F del puerto con la isla Libertad) empleando 100 Ci de THO como traza- 
dor. El trazador se diluyó en 200 litros de agua de la bahía. Se inyectó en forma continua, 
mediante una bomba peristáltica, durante 5 horas, en el rastro de la barca utilizada. 


-Marcado del vertido de aguas servidas del colector 12 de diciembre con '*'L 


En la mayor parte de los trabajos con radioisótopos que se describen en este artículo par- 
ticiparon expertos de la sección de hidrología isotópica del OIEA, quienes contaron con 
sus propios dosimetros. Esto permitió comparar las dosis informadas por el OIEA con las 
dosis determinadas localmente, en la Dirección Nacional de Tecnología Nuclear. 

Para cada trabajo de campo se solicitó y se obtuvo, de parte de la autoridad competente, la 
correspondiente autorización para el uso de radioisótopos sobre la base de la justificación 
de la práctica, comparada con posibles alternativas, su optimización y la limitación de la 
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dosis tanto para los trabajadores ocupacionalmente expuestos como para el público en 
general. Durante las operaciones se aplicaron medidas de seguridad y se realizó un moni- 
toreo ambiental, por los oficiales de protección radiológica, de acuerdo con las normas. 
Se establecieron medidas de seguridad para enfrentar situaciones no planificadas que 
involucraran riesgo de exposición o de contaminación. Se produjo un incidente durante 
las maniobras previas al marcado de la pluma de descarga del emisario de Punta Brava 
con **!T, que involucró a un funcionario del OIEA. Se resolvió satisfactoriamente aplican- 


do las medidas de seguridad previstas para el caso. El resto de los trabajos se desarrolla- 
ron sin incidentes ni accidentes. 


Para el procesamiento de los datos y su representación matemática, la bahía se subdividió 
en celdas a cuyos puntos medios (correspondientes a estaciones de muestreo) se les asig- 
naron perfiles de temperatura, conductividad y concentraciones de tritio (Figura 3). 

Además, se midieron las concentraciones de tritio en puntos de muestreo distribuidos a lo 


largo del cinturón de playas situadas al Este de la boca de la bahía, hasta la desembocadu- 
ra del arroyo Carrasco. 
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Figura 3. Sub-división de la bahía en celdas. 
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La Figura 4 presenta algunos detalles batimétricos uilizados para programar las 
operaciones e interpretar los resultados obtenidos en los trabajos de campo, expresados en 
un formato adecuado para implementar una simulación digital del campo de velocidades 
y de niveles del agua variables asociados a las mareas astronómicas y meteorológicas. 
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Figura 4. Detalles batimétricos relevantes para los trabajos de campo y subdivisión del dominio 
espacial para los trabajos de simulación digital (CFD) mediante aplicación del método de elementos 
finitos. Un ejemplo de cálculo numérico aparece en la Figura 5. 


Figura 5. Algunos resultados de la simulación digital del campo de velocidades durante la marea 
creciente obtenidos por el IMFIA, Facultad de Ingeniería, Universidad de la República. 
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Para investigar un posible efecto de contaminación de las playas debido a la liberación de 
contaminantes de la bahía en condiciones de marea descendente (o sea con carreras de 
marea hacia el Este) se inyectaron en forma súbita 100 Ci de THO en un lugar que apare- 
ce indicado en la Figura 6: cerca del extremo distal (respecto de la costa) de la escollera 
Sarandí, en las proximidades y al Este del canal de acceso al puerto de Montevideo. 


Durante las campañas de marcado y estudio de la distribución de las aguas contaminadas 
que escapan de la bahía de Montevideo se efectuó un estudio sobre cinco transectas fijas. 
En cada transecta se definieron seis puntos de muestreo, separados por 1 km, y comen- 
zando siempre con un punto en la costa [3], [4]. La Figura 5 muestra las transectas men- 
cionadas, donde se determinaron concentraciones de tritio, conductividad y temperatura 
durante tres semanas. 
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Í 0,09610.004 ¡CUm' 
| O 1/4/93 a las 9:10h 
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5 24/93 a las 16:55h | 
0,05540.003 Cin - ; ! 
28/3/93 a las 16:05h 0IAA10.003 paCas Figura IL1, Resultados de los análisis de tritio de muestras medias | 
350 0.05610.003 CUA? 28/0/93 a las 15:35h _ representativas de la totalidad de cada transecta medidas con enriquecimiento 
y 1/4/93 u las 11:300 | electrolítico previo. | 


Figura 6. Ubicación del punto donde se inyectaron en forma súbita 100 Ci de THO y ubicación de 
las transectas donde se determinaron concentraciones de tritio, junto con perfiles de conductividad y 
temperatura durante tres semanas. 


Durante las campañas de marcado y estudio de la pluma de descarga del emisor subacuá- 
tico de Punta Brava se midieron las actividades del trazador en puntos de la pluma de 
descarga del emisor, que se movía arrastrada por las corrientes [3]. En la referencia [9] y 
en su bibliografía se pueden encontrar detalles sobre las técnicas empleadas. 


Se efectuaron dos campañas. En cada una se inyectó una solución de ioduro de potasio 
marcada con **!I en una de las piletas de tratamiento previas al emisor, durante 2 h 45 
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min, a razón de 0.5 l/min. La concentración del trazador era 1377 u-Ci/m'. Se siguió el 
despliegue de la pluma de descarga marcada durante ocho horas cada vez, a tiempo real, 
utilizando detectores de centelleo. La pluma fue detectada y medida en cinco ocasiones, 
hasta que la dilución del trazador impidió continuar con las mediciones. Se utilizaron 
biplanos y uranina para posicionamiento y mediciones de velocidad del agua. 


Durante las campañas de simulación de un accidente de contaminación en la playa Poci- 
tos, se inyectaron en forma localizada y súbita soluciones de uranina y se midió el trans- 
porte advectivo y la dispersión de las manchas de colorante. Se repitieron los experimen- 
tos en las cercanías del emisor subacuático. Al principio, se aplicaron técnicas de medi- 
ción a bordo de las concentraciones de trazador, en condiciones de flujo continuo median- 
te un fluorímetro, durante una secuencia de maniobras de navegación efectuadas sobre la 
mancha en cada una de las etapas de evolución seleccionadas. Estas técnicas no permiten 
suficiente resolución espacial para estudiar el mencionado efecto de histéresis y detalles 
debidos a los efectos de micro-circulación. Por este motivo se decidió registrar la mancha 
de trazador, por medios ópticos y desde el aire, para estudiar posteriormente la estructura 
de sus bordes mediante herramientas de análisis de patrones bidimensionales. Con el 
propósito de disminuir lo más posible la histéresis, se diseñó y construyó un corral para la 
inyección de trazadores [5], que resultó muy eficaz, por lo cual fue utilizado en las cam- 
pañas siguientes (Figura 7). 


Figura 7. Corral concebido y construido por Alejandro Malek para inyectar trazadores de modo de 
minimizar efectos de histéresis y facilitar la investigación de los efectos de escala en aguas costeras. 


Para determinar los parámetros de la dispersión se efectuaron dos inyecciones en la capa 
superior del cuerpo de agua estratificado y una inyección en el cuerpo no estratificado. Se 
describió el efecto de escala suponiendo un coeficiente de dispersión horizontal del tipo 
E, = A-r* donde q representa el parámetro de efecto de escala. Para estimar q se aplicó 
el método de los circulos equivalentes [5]. Para estimar la cantidad de trazador a inyectar 


en la capa superior del cuerpo de agua se asumió a priori q=1. Estos aspectos se retoman 
en la segunda parte del artículo. 


Durante las campañas de dinámica de sedimentos se efectuó el seguimiento del material 
dragado marcado y descargado por la draga. 


Revista de la Soc. Científica del Paraguay, 22 (1):13-46 ISSN 0379-9123 


Se determinaron parámetros de dispersión de la nube de partículas en suspensión, se esti- 
maron tasas de sedimentación y se estudió el movimiento del material sedimentado. 

Se efectuaron mediciones relacionadas con la profundidad náutica y con la optimización 
del tiempo de bombeo en dragas de succión [7]. 


Durante las cuatro campañas destinadas a estudiar la capacidad de autodepuración de las 
aguas costeras de Montevideo en la región de Punta Lobos, en todos los casos se inyectó, 
en la zona prevista para el emisor subacuático, en forma súbita y puntual, una masa de 
efluente conteniendo la población bacteriana, acompañada de un trazador fluorescente y 
conservativo (rodamina WT) para determinar la dilución por el proceso de advección- 
dispersión. 

El objetivo final fue obtener una estimación del 7, bacteriano en diferentes escenarios. 


Se utilizó uranina como trazador secundario cualitativo, para visualizar la ubicación del 
vertimiento efectuando marcados sucesivos al finalizar cada etapa de mediciones. En 
cada etapa se tomaron muestras para determinar unidades formadoras de colonias y con- 
centración del trazador en laboratorio. 

Además se determinó la posición de cada punto de toma de muestras, los perfiles vertica- 
les de velocidad de corriente, de temperatura y de conductividad, la velocidad local del 
viento, la turbidez del cuerpo de agua, el pH, la DBO y otras variables características de 
cada escenario en el que se produce el transporte advectivo-dispersivo y la inactivación 
bacteriana. 


Resumen de resultados experimentales 


La dinámica del agua en la bahía se puede simplificar asumiendo la presencia de dos 
flujos separados, cuyas características dependen de la ubicación variable del extremo de 
la cuña salina en las aguas costeras de Montevideo: un flujo de agua de baja conductivi- 
dad que ocupa una capa de menos de 4 m de profundidad y un flujo de agua de conducti- 
vidad más elevada que ocupa las zonas con profundidades mayores a los 4 m. 

Este último flujo penetra a través del canal de acceso al puerto. El otro flujo penetra pre- 
ferentemente a través de la boca Oeste. 

La estratificación del agua en la bahía es permanente, lo cual debe ser tomado en cuenta 
al construir modelos matemáticos para la simulación de la dinámica y la renovación del 
agua. 

Durante la marea creciente (carrera de marea dirigida hacia el Oeste) el patrón de flujos 
en la bahía es variable, dependiendo de los vientos, pero en promedio corresponde a una 
fase de llenado. 

Durante la marea decreciente (carrera de marea dirigida hacia el Este) el flujo que des- 
ciende por el estuario, próximo a la costa, es desviado por la escollera Oeste y se divide 
en dos. 

Se produce entonces un flujo secundario que pone en movimiento a una parte significati- 
va del agua de la bahía bajo la forma de un flujo circular en el sentido de las agujas del 
reloj: este movimiento es el principal responsable del mezclado y de la renovación del 
agua (Figura 8, izquierda). 
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periodos de marea decreciente. 

Figura 8. Izquierda: flujo circular en la bahía, cuando la carrera de marea se dirige hacia el Este. 
Derecha: estimación de un tiempo de renovación del agua de la bahía a partir de los datos de activi- 
dad total remanente (en coordenadas logarítmicas) versus tiempo. En este caso el resultado fue de 
29 horas. 


La marea decreciente corresponde a una fase de vaciado: la bahía se comporta como una 
bomba de descarga intermitente de contaminantes en las aguas costeras adyacentes. 

A partir del volumen y la concentración en cada una de las celdas en las que se dividió la 
bahía (Figura 3) se puede estimar la actividad de tritio remanente en función del tiempo 
(Figura 8, derecha). 

Considerada en conjunto, el agua de la bahía se renueva con un tiempo promedio de 29 
horas, si se asume un único compartimento representativo. Si se consideran los compar- 
timentos parciales definidos por las celdas de muestreo, el tiempo de renovación disminu- 
ye al pasar de las celdas de la parte suroeste hacia las celdas de la parte noreste. En la 
mayor parte de esas celdas el tiempo de renovación es relativamente corto. Las zonas de 
renovación más lenta corresponden a la desembocadura de los arroyos Miguelete y Pan- 
tanoso, a las dársenas del puerto y a la región comprendida entre el espigón F y la central 
termoeléctrica. 


En las campañas efectuadas en marzo y abril de 1993 se encontró un cuerpo de aguas 
costeras con una estratificación bastante neta. Generalmente se observaron dos capas de 
agua de diferente salinidad, separadas por una picnóclina: una capa superior, de agua más 
dulce, de río, y una capa inferior, de agua más salada, proveniente del océano. 

Al mismo tiempo la conductividad variaba horizontalmente aumentando hacia el Este 
tanto en la superficie como en el fondo y la profundidad de la picnóclina aumentaba hacia 
el Oeste. 

Una diferencia de conductividad entre el fondo y la superficie de 10 mS/cm a través de 
una columna de agua de 6 m, con una capa de 1 m de espesor a 3 m de profundidad, don- 
de la conductividad varía bruscamente en 5 mS/cm (haloclina), es representativa de la 
situación. Las diferencias de temperatura medidas (menos de 2 * C) entre la superficie y el 
fondo del estuario no inciden en forma significativa en la estratificación, puesto que 2*C 
corresponden a una diferencia de densidad de 0.6 g/l, contra los 7 g/l correspondientes a 
una diferencia de conductividad de 10 mS /cm. 


Revista de la Soc. Científica del Paraguay, 22 (1):13-46 ISSN 0379-9123 


De un día a otro se producían variaciones en más o en menos respecto de los valores de 
conductividad obtenidos en el mismo punto y a la misma distancia de la superficie. Estas 
modificaciones, bien conocidas, se deben al avance o al retroceso de masa de agua salada 
o de agua dulce. Pero el patrón de estratificación vertical se mantuvo en todos los casos. 
Durante los días 24 al 28 de mayo se produjo una invasión de agua de baja salinidad pro- 
veniente del Oeste. La conductividad en la playa Pocitos, durante esos días, no descendió 
a los valores bajos (4 mS/cm) medidos en la playa Ramírez. A partir del día 28 se produjo 
una invasión de agua más salina (conductividades superiores a 35 mS/cm) proveniente del 
Este. 


Las estimaciones de velocidad efectuadas utilizando biplanos y las mediciones con co- 
rrentómetros muestran en general una diferencia significativa en las velocidades prome- 
dio entre la capa de agua predominantemente dulce y la capa de agua predominantemente 
salada. Esas diferencias dependen del instante en el ciclo de marea astronómica y de los 
vientos predominantes. 

En la zona del emisor subacuático, durante el segundo experimento con !%'I, las velocida- 
des en la capa superior fueron aumentando desde 140 m/h hasta 310 m/h, mientras que las 
velocidades en la capa inferior aumentaron desde 360 m/h hasta 1200 m/h, para descender 
luego a 960 m/h hacia el final del experimento. 

El agua circula a lo largo de una franja cuyo ancho varía entre 1 y 5 km, con la región 
más estrecha a la altura de Punta Carretas. Los procesos de transporte son muy rápidos a 
lo largo de la costa: en dos o tres días una partícula que parte próxima a la boca de la 
bahía puede alcanzar la playa Carrasco, a más de 15 km de distancia. 


La primera detección de tritio se efectuó a las 2 horas de ocurrida la pleamar siguiente a la 
inyección del trazador. En ese entonces el agua marcada había llegado ya a la transecta 1 
(playa Ramírez, ver Figura 6) con una concentración elevada en un punto próximo a la 
costa. 

Al día siguiente no se encuentra tritio en la transecta 2 (playa Pocitos) ni en la zona más 
próxima a la costa de la transecta 3. Pero se halla el trazador en la zona más alejada de la 
transecta 3 y en la transecta 4. Esto indica que en la playa Pocitos se produce un efecto de 
trampa que tiene como consecuencia que las descargas locales de contaminantes adquie- 
ran más importancia allí que en otras zonas de la costa. 


La aparición de agua de conductividades muy diferentes, marchando en sentidos opuestos 
y mezclándose con el agua marcada no llegan a producir una renovación apreciable del 
agua de la zona costera investigada: el tritio continuó estando distribuido a lo largo de la 
cadena de playas, a concentraciones más bajas pero medibles sin dificultad. 


Durante los estudios de estabilidad de la pluma de descarga del emisor de Punta Brava, la 
sensibilidad del detector en las condiciones del experimento fue de 2.1 cuentas por se- 
gundo por metro cúbico y micro-curio. Efectuando un balance de masa se puede estimar 
una concentración de 442 u-Ci/m! a la salida del emisario, para un flujo de 1 metro cúbi- 
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co por segundo. Esto corresponde a una respuesta de 998 cuentas por segundo en el detec- 
tor. 

Como el máximo valor medido en el penacho formado a la salida de los difusores fue de 
181 cuentas por segundo, esto indica que el factor de dilución en ese punto era de 5.1. 
Este coeficiente de dilución inicial determinado con trazadores es casi la tercera parte del 
coeficiente estimado mediante la teoría de chorros en medios estratificados [3]. 


La medición de la actividad del trazador mostró que el agua descargada quedaba atrapada 
por debajo de la picnóclina, en forma completamente estable. Viajaba hacia el Oeste, 
siguiendo la carrera de marea ascendente, aumentando paulatinamente su profundidad y 
expandiéndose horizontalmente más y más. 

Cuando las velocidades de flujo son bajas, las aguas residuales viajan dentro de una capa 
de menos de 1 m de espesor que se ensancha progresivamente a medida que el vertido se 
aleja de los difusores. 

A velocidades mayores aumenta el transporte por difusión turbulenta vertical en el fondo 
del estuario, por lo que se produce un mezclado vertical del agua situada entre el fondo y 
picnóclina. Se detecta experimentalmente a través del perfil vertical de actividad del tra- 
zador en la 5* detección de la pluma (Figura 9). 


Tiempo: 11:44 a 12:06 h 
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Figura 9. Izquierda: segunda detección de la pluma. Derecha: quinta detección de la pluma. Las 
líneas de iso-actividad para la quinta detección de la pluma están corregidas para tener en cuenta 
efecto del transporte advectivo durante el intervalo de tiempo involucrado en las mediciones. 


La dilución del trazador aumentó muy rápidamente y las concentraciones alcanzaron en 
algunos puntos valores por debajo de los límites de detección. No obstante, la halóclina 
no se desvanece: permanece como un techo difusivo para las aguas residuales vertidas por 
el emisario, debido a lo cual la transferencia de agua contaminada hacia las capas superio- 
res se produce con dificultad. 


El espesor inicial de la pluma de descarga del emisor sub-acuático puede estimarse a 
partir del factor de dilución minimo, del flujo descargado por unidad de longitud del tra- 
mo de difusores y de la velocidad de la corriente transversal al tramo de los difusores. Si 
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se toma una velocidad de 400 m/h (que corresponde al comienzo del experimento), una 
longitud de 67 m para el tramo de difusores y el flujo total descargado, se obtiene un 
espesor de 0.7 m para el inicio de la pluma [3]. 

El espesor estimado a partir de los perfiles de actividad fue de 2 m. 

Estas dos estimaciones son compatibles por dos motivos: (1) el trazador no viaja en una 
lámina de concentración constante y (2) el detector es capaz de captar fotones gama pro- 
veniente de puntos relativamente alejados. 

Aunque el peso del aporte de un elemento de volumen al número total de cuentas cae 
rápidamente con la distancia entre dicho elemento y el detector, este último comienza a 
producir una señal significativamente por encima del ruido de fondo cuando todavía se 
encuentra algo alejado de la nube de trazador. 


En las campañas llevadas a cabo en la segunda mitad del mes de mayo y en el mes de 
junio de 1995, la estratificación vertical típica de la época estival había desaparecido, y la 
diferencia de temperatura entre la superficie y el fondo, de dos o tres de grados se había 
reducido a menos de un grado Celsius. 

Se obtuvieron valores del parámetro de efecto de escala próximos a 1 al comienzo del 
proceso. Los valores del parámetro de escala fueron disminuyendo continuamente hasta 
anularse cuando la nube de trazador se había hecho lo bastante extensa, aproximadamente 
luego de transcurrida una hora y media desde el momento de la inyección. 


Cuando la mancha es lo suficientemente grande (su menor dimensión superior a los 
700m) se puede utilizar el modelo Gaussiano para una caracterización paramétrica del 
proceso de dispersión del trazador. Para los coeficientes de dispersión longitudinal se 


A on 2 y 
obtuvieron valores próximos a 1] mi ( K, 20.83 M 2 ) y cerca de veinte veces meno- 
Ss Ss 


res para los coeficientes de dispersión transversal (K. = 0.018 m/ ). 
S 


En las campañas de 2001 destinadas a estudiar la capacidad de autodepuración de las 
aguas costeras en la región de Punta Lobos se estimaron, para el parámetro equivalente al 


Tp , valores comprendidos entre 6 y 12 horas, tomando valores mínimos en los escena- 


rios correspondientes a las dos últimas campañas y valores máximos en los escenarios 
correspondientes a las dos primeras. Si bien las temperaturas del aire durante las últimas 
dos campañas eran relativamente elevadas para la época del año, la temperatura del agua 
no lo reflejaba, como cabía esperar. En las primeras dos campañas la turbidez del cuerpo 
de agua era significativamente mayor que la turbidez durante los dos últimos trabajos de 
campo. 


SEGUNDA PARTE: MODELOS MATEMÁTICOS 


En general conviene distinguir dos tipos de modelos matemáticos: los que se utilizan para 
el diagnóstico y los que se utilizan para el pronóstico de las propiedades de un cuerpo de 
agua. Ambos tipos de modelos matemáticos, de diagnóstico y de pronóstico, se pueden 
combinar con el empleo de trazadores. 
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Los modelos de diagnóstico permiten caracterizar los procesos de transporte de masa 
mediante un conjunto de parámetros significativos. 

Los modelos de pronóstico permiten la predicción de algunas características de importan- 
cia de los procesos de transporte, basada en los parámetros hallados aplicando los 
modelos de diagnóstico y en otra información relevante. 

Algunas veces un mismo modelo se puede utilizar tanto para diagnóstico como para pro- 
nóstico. 

Actualmente, con fines de pronóstico, la mayor parte de las veces se utilizan sofisticados 
códigos de cálculo en el marco de la CFD, que se calibran para cada aplicación específica. 
Muchos de estos códigos están disponibles bajo la forma de paquetes de software comer- 
ciales. Permiten obtener soluciones numéricas válidas para determinado conjunto de valo- 
res de los parámetros. Los procedimientos para graficar las soluciones, como el campo de 
velocidades que muestra la Figura 5, ayudan a identificar patrones interesantes como los 
vórtices horizontales en la región nordeste de la bahía que muestra esa figura. 

Ahora bien, no siempre se pueden hallar relaciones de validez general o se pueden identi- 
ficar procesos interesantes a partir de un conjunto de simulaciones digitales si no se dis- 
pone de información adicional para guiar el cálculo numérico. Si resulta factible (cosa 
que muchas veces no lo es) el recurrir a modelos analíticos permite obtener resultados 
para todos los valores de los parámetros, facilita la extracción de relaciones de validez 
general y la comprensión de esas relaciones, y puede suministrar información útil para la 
simulación digital. 

Los modelos matemáticos que se describen y discuten a continuación son algunos de los 
modelos analíticos utilizados para diseñar e interpretar los experimentos con trazadores 
reseñados en la primera parte del artículo. 


Modelos de compartimentos para describir la renovación del agua en la 
Bahía de Montevideo 


En principio se podría plantear un modelo matemático de compartimentos múltiples para 
la cinética del tritio en la bahía. 
Supongamos que la bahía se subdivide en p compartimentos, que se supone son homogé- 
neos para los fines del modelo. 
Las ecuaciones de balance de los volúmenes de los compartimentos se pueden escribir 


: d p P , 
asi: == 0; A Oj +00 ¿+04 1=1L2,..¿P (1) 


dt jah jkió  j=L ¡ei 
Aquí O, ¡es el influjo de volumen que recibe la celda i desde la celda ¡, V, es el volu- 


men de la celda ¡, O 


ci 


es el eflujo de volumen desde la celda ¿ hacia las aguas costeras, 
Q,, es el influjo de volumen hacia la celda ¡ desde las aguas costeras y O es la suma de 


todos los influjos de agua de arroyos o vertidos sanitarios que van a parar a la celda 7. 
En general los flujos hacia y desde las aguas costeras son nulos, excepto en las caldas 
adyacentes a esas aguas. Los influjos de agua dulce son diferentes de cero en las celdas 
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que reciben vertidos sanitarios o son adyacentes a la desembocadura de cursos de agua 
como los arroyos Miguelete y Pantanoso. 

Las variaciones de volumen en las condiciones en las que se llevaron a cabo los trabajos 
de campo no superaron el 10 % del volumen de cada celda. Por ello y con el fin de sim- 
plificar el modelo, supondremos que los volúmenes permanecen constantes. Una discu- 
sión sobre este punto puede hallarse en la referencia [10]. 


Supongamos ahora que la actividad de tritio en la celda / se representa por A (e). Un 
balance de actividades en una escala de tiempo en la cual el trazador radioactivo se puede 


considerar conservativo (del orden de cientos de horas comparada con los más de 10 años 
de tiempo de semi-desintegración del tritio) conduce al siguiente sistema de ecuaciones 


diferenciales ordinarias en términos de las probabilidades pl) 40 de hallar una 


partícula de trazador en cada una de las celdas en el instante ;: 


Cp l)=Em PO i=L2p 0) 
dt da 


P 
20; +0 
daniel 


v. 


J i 
Se supone que no hay reflujo de THO desde las aguas costeras hacia la bahía. Esta supo- 
sición se basa en los resultados obtenidos en los trabajos de campo. 


EH) MB) e) 


Por definición m,; = =2 cuando ¡4 ¡. Cuando ¡= j, m,;= 


En forma matricial: 

Aquí la matriz cuadrada M = Im : sl y P(t) es un vector columna de px1 formado 
PXp 

por las probabilidades P (e) = 44e) La solución Ple) de la ecuación matricial (3) para 


cualquier condición inicial P(0) se puede construir combinando linealmente soluciones 


Aeotos e Ñ > a 
de la forma e”* -v, donde y) ¿es uno de los valores propios y V, es un vector propio 


correspondiente a ese valor propio. 
Siempre se puede asumir que los valores propios son diferentes dos a dos sin modificar la 
bondad del ajuste a los resultados experimentales [10]. 


En principio, esta estructura de modelo se podría aplicar con p=20 compartimentos (un 


compartimento por cada celda), pero los datos experimentales disponibles no permiten la 
construcción de un modelo robusto con 20 variables de estado y centenares de paráme- 
tros: son demasiado escasos y ruidosos, no permitiendo una adecuada resolución temporal 
y espacial, aún si se añadiera información proveniente del modelado por métodos de CFD 
del tipo utilizado en paralelo con los trabajos de campo. 
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En lo que sigue se consideran un modelo clásico con p=2 y un modelo con p=3. Este 
último no se reduce al modelo general de p compartimentos homogéneos porque el ter- 


cer compartimento es fractal. 

No obstante, con dos o tres compartimentos no se puede dar cuenta de ciertos fenómenos 
observados durante los trabajos de campo. En particular, la evolución temporal de las 
actividades en algunas de las celdas ubicadas en el Norte de la bahía exhibe un patrón 
característico, con incrementos y decrementos sucesivos y cada vez menores. Esto parece 
estar relacionado con la marcada asimetría en los flujos de intercambio entre celdas veci- 
nas, asociado con el flujo circular del agua durante la marea decreciente que se ve en la 
Figura 8 (izquierda). 

En general el fenómeno mencionado se manifiesta a nivel de la estructura del modelo 
matemático a través de la aparición de pares de valores propios de la matriz 


M= Im, al que son complejo-conjugados con parte real negativa y parte imaginaria no 
Jlpxp 


nula. 

La recta de mejor ajuste que aparece en la Figura 8 (derecha) se trazó excluyendo los 
datos provenientes de las mediciones efectuadas durante las primeras 18 horas a partir de 
la inyección del trazador. Pero cuando la recta se extrapola al instante de la inyección, se 
obtiene una actividad inicial de 133,3 Ci, superior a los 100 Ci inyectados. 

Este hecho junto con los detalles espaciales y temporales de la cinética del tritio en la 
bahía, los datos batimétricos disponibles, información sobre el campo de velocidades del 
agua en el interior de la bahía obtenida con biplanos, correntómetros, trazadores y resul- 
tados de simulación con métodos de CFD (modelos numéricos a elementos finitos pro- 
mediados en profundidad), sugieren que la bahía puede ser dividida en dos grandes com- 
partimentos: Nord-Este (NE, con un volumen promedio de 1,6x10” m*) y Sud-Oeste (SW, 
con un volumen promedio de 3,2x107 m?). 

Teniendo esto en cuenta, se construyó un modelo clásico de dos compartimentos, NE y 
SW, con las actividades de tritio en cada uno de esos compartimentos como variables de 
estado. 

La Figura 10 muestra los dos compartimentos, los flujos de volumen que los interconec- 
tan entre sí y con las aguas costeras, y los aportes de agua dulce proveniente de tierra. 


l Qs IE 


Q 
Qe (SW) | Q | (NE) 


Figura 10. Esquema de un modelo de dos compartimentos para la cinética del tritio. 


Si A (2) y Á (1) son las actividades en los compartimentos SW y NE respectivamente, el 
balance de actividades se reduce a las dos ecuaciones siguientes: 
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Ae (O.+01) , Q 


= Á 4 
a v, E (4 a) 
NS Yara, (4b) 
dt v, V 
Con la condición inicial A(0) =0 y A2(0)= Ap la solución es la siguiente: 
0, 1 
Aj(t)= 4: : «[exp|-|A ¡|-t]-exp|-(2,|-1 (5 a) 
(4) 0 V, (ea) ( [ | 1] ] [ | Al ) 
Ar(i)= 4, (+) (e 2| 2) expl-|2,]-e] (a 2) exp-|2 4]: 2) 59) 
Aquí 4, =-—|2 | and 4,=-4)| (con [2,|>|2 ,|) son las raíces siempre negativas de la 
ecuación característica: q (e +2.) + da A+ La O =0 (6) 
Y, e E 


Sumando miembro a miembro (5 a) y (5 b) se obtiene la actividad total en la bahía en 
función del tiempo a partir de la inyección del trazador (que se supone instantánea): 


(7) 


(2-21) 


Alt)= 4 na expl- Ay]: e|- l2,]- expl- 22) a) 


La actividad inicial A, corresponde a los 100 Ci distribuidos en el compartimento nores- 


te, que se comunica directamente solo con el compartimento suroeste. 

La recta que se dibujó a la derecha en la Figura 8 se ajusta mejor a los datos que corres- 
ponden a las últimas etapas del vaciado de tritio de la bahía pero no es representativa de la 
etapa Inicial. 

Por el contrario la ecuación (7) describe bien esta primera etapa del proceso. 

Durante esta primera etapa el tritio ocupa fundamentalmente el compartimento NE y va 
siendo transportado hacia el compartimento SW. El compartimento SW comprende la 
boca de la bahía y se comunica directamente con las aguas costeras: el tritio solo puede 
abandonar la bahía a través de este compartimento. Por este motivo la actividad 4 Al) en 


el compartimento NE decrece a partir de 4, mientras que la actividad Ale) en el com- 


partimento SW crece desde 0, acompañada de eflujo de THO hacia las aguas costeras que 
es pequeño en relación con el influjo de THO proveniente del compartimento NE. 
El ajuste de parámetros, que se describe en la referencia [10], arroja 0,117 hr"! para |2,| y 


0,034 h” para |A,|. 
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Una vez desaparecido el efecto de 2 ,|-exp[- 2 2]-:), la actividad total remanente se 


puede describir por la recta asintótica: log, A(t)= la ]-1+log, EN (8) 
y) 
1 - 1 
4) 
Esta expresión corresponde a la recta de la Figura 8. Entonces: Az (0) = a. (9) 


e) 


Como por (9) Az (0) es mayor que la concentración inicial de tritio 4, , esta última queda 


por debajo del punto de intersección de la recta con el eje de los tiempos. 
La actividad total A(t) viene representada por curva de la Figura 11. 


A PP E 
O. 0 A AI, CE 
POEMA: TERR CRA is inte: 
ls o 


á 


: 
| 
| 
| 


Figura 11. En abscisas el tiempo en horas, a partir de la inyección del trazador. En ordenadas, en 
escala logarítmica, el porcentaje de actividad de tritio. Los puntos corresponden a mediciones. La 
curva ajustada corresponde a la actividad total de acuerdo con la fórmula (7). 


2 
Al inicio A A(0)= 0 y = log, A(t)< 0. Luego, hasta unas 70 u 80 horas a partir de la 
É É 
inyección, log, Alt) se comporta aproximadamente como lo indica la fórmula (8). 
Si Í es una variable aleatoria que representa el tiempo que reside una partícula de traza- 
dor en la bahía, la probabilidad PIT >+1| se puede estimar recurriendo a la denominada 


función de remanencia R(r) = = del trazador en la bahía: Pl? > 1|= R(+) 
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La función de distribución f; (+) de la variable aleatoria 7 se puede hallar así, teniendo en 


cuenta que PIT <t|=1-P|F >+1|: ro rr <: 2 la (0) 
E dt Ay dt 
+00 A(t) 


Teniendo a su vez en cuenta que e fo (0) di= 6 di el valor promedio 


L= l E a (1). dt del tiempo de permanencia (o tiempo de renovación) del trazador se 


rg (11) 


Para el modelo clásico de dos compartimentos obtenemos: 


puede estimar así: 


pd 
Pal 14 
Sustituyendo en (12) los valores 0,117 h”' para [2,| y 0,034 h para |2,| resulta 


(12) 


T, =37h , mayor que las 29 horas obtenidas a partir del modelo de un compartimento. 


Si bien el ajuste de la formula (7) a los datos experimentales es bueno durante la mayor 
parte del proceso de renovación del agua de la bahía, el curso temporal tanto de la cola de 
si Alt . Ai Az 
la función R(+)= al) como de las colas de las funciones 410) y 440) 
Ao Ao Á 
actividades en los compartimentos y la actividad total se desvanecen más lentamente que 
una función exponencial del tiempo. 
Esto plantea un problema que no puede resolverse en el marco del modelo lineal de dos 
compartimentos clásicos. Si se mantuviera la validez de este modelo, asintóticamente para 
dde ; , A 
t > +00 la actividad en la bahía debería desvanecerse como 4) - a i s | «exp[-|2| | 
21141 
pero no es esto lo que parece que ocurre. El lento desvanecimiento de la actividad apunta 
a un efecto de dispersión anómala de las moléculas de THO, posiblemente vinculado a 
caminatas aleatorias que se encuentran severamente restringidas en algunas partes del 
sistema. 


sugiere que las 


La distribución de zonas muertas en el interior de la bahía junto con la geometría del 
fondo, incluyendo la distribución del barro fluido, sugiere una posible solución al pro- 
blema: introducir, además de los compartimentos clásicos NE y SW un tercer comparti- 
mento, fractal. Con esto abandonamos el marco de la teoría clásica de compartimentos y 
debemos recurrir a la teoría no convencional. 

Actualmente hay dos enfoques para la teoría no convencional de la cinética de comparti- 
mentos: el fraccionario [11] y el fractal [12]. 

El enfoque del análisis de compartimentos fraccionario y lineal se basa en introducir flu- 
jos hacia o desde un compartimento fractal proporcionales a derivadas de orden fraccio- 
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nario de las concentraciones. Eso se hace en la referencia [13] para la cinética del tritio en 
la Bahía de Montevideo. 

El enfoque alternativo del análisis de compartimentos lineal no convencional, denomina- 
da análisis fractal, se basa en introducir coeficientes cinéticos que son funciones decre- 
cientes del tiempo ke), y tienden a cero para f —>-+00 como funciones potenciales del 


tipo £7P para O< p<l. 

Ambos enfoques permiten tomar en cuenta procesos de transporte restringidos en medios 
heterogéneos como el compartimento fractal de la bahía, pero no son completamente 
equivalentes. 


En la referencia [14] se construye un modelo de tres compartimentos no convencional que 
generaliza a y es compatible con ambos enfoques: fraccionario y fractal. La Figura 12 
presenta un esquema del nuevo modelo: 


Figura 12. Modelo de tres compatimentos, uno de ellos fractal (el compartimento 3). 


Utilizando el enfoque fractal y aplicando un método de perturbaciones se obtiene la si- 
guiente solución aproximada, válida para la última etapa del proceso de desaparición del 
trazador [14]: 


Lope 
lp ia lo po sd Kkle Solo) de 
(1-7, [0 da) sel Hu)e” alo) (13 a) 
l pro -Az(0) 
l+8, (Es ku) e" <a) 


a0)> 
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En la ecuación (13 a) 4 ee) es la suma de las actividades de tritio en los compartimentos 


NE y S que forman un compartimento combinado R. Los volúmenes de los comparti- 
mentos fractal F y NE+SW combinados se Eo por Y, y V, respectivamente. El 


parámetro positivo sin dimensiones «, == 
FOR 


se asume es pequeño respecto de 1. 


En la ecuación (13 b) 4 E) es la actividad de tritio en el compartimento fractal F. En 
ambas ecuaciones (13), 4g(0) viene dada por la fórmula (9). Por definición: Az =|2 | 

La actividad total de tritio se obtiene sumando ambas ecuaciones miembro a miembro. 

En la primera etapa del proceso de desaparición del trazador (que se produce al unísono 
en ambos compartimentos homogéneos, NE y SW) 4 ale) viene disminuyendo y es mayor 
que Á she) en al menos un orden de magnitud. A su vez A ple) viene aumentando a partir 


de cero. 
El flujo de THO se dirige desde R hacia F. Esto se representa por el término proporcional 


a (E [0 e «en (13 a). La actividad en R disminuye más rápido que 
R 


ett porque al eflujo de THO hacia las aguas costeras (supuestamente sin reflujo) se 
añade un eflujo neto de THO hacia el compartimento fractal. 


El factor es, f k(u).du termina por anularse y luego se hace negativo. Entonces 
V, 


1 j : k(u). du) e7'«** pasa por un valor mínimo negativo, que cabe esperar sea muy 
VD, 


pequeño en valor absoluto respecto de 4, (0), y luego tiende a cero por valores negativos 

muy pequeños respecto de los valores que toma durante la etapa inicial. El otro término, 

siempre positivo, que aparece en el numerador de (13 a), pasa a ser dominante. Se invierte 

el sentido del flujo entre F y R, que ahora se dirige desde F hacia R. En esta última eta- 
Í ÁLt 

pa Ez es mayor que —* 

F R 
mantienen del mismo orden de magnitud durante el resto del proceso. 


, pero esas concentraciones promedio de actividad se 


Cuando £-—>-+00 la influencia de (E J ' Hadas) ee sobre la evolución de 
Va 70 


A glt)se hace despreciable. Ambas actividades, A p(t) y A pt) se desvanecen al uní- 
sono, al igual que la actividad total Alt) , así: 

4 lito du d Aito du 
) e =—€ 


kt 
dt 


(14) 
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De la fórmula (10) para la función de densidad de probabilidad f (+) para el tiempo 


aleatorio de residencia de una partícula de THO en la bahía y de la fórmula (14) se des- 


4) feto du 
prende que para Í >+00: folt)o eE (15) 


Si el coeficiente cinético que caracteriza el transporte hacia y desde el compartimento 
-p 
fractal (F) sigue una ley potencial k(t)=k, ( + 2) (0< p<1) se desprende la siguien- 


te relación asintótica durante la última etapa de la desaparición del trazador: 


log, A ¿(t) = (const) puo5 7) AUS 0<p<l (16) 


api Ain 


En este caso la densidad de probabilidad de residencia tiende asintóticamente a cero como 


ON 
una exponencial dilatada: ¿(t)oc e Y AS O<p<l a? 
Í Pp 


En la Figura 13 se resumen todas las etapas del proceso de desaparición de THO en una 
curva ajustada a los datos experimentales disponibles [14]. 


2 
> 
- 
ha 
8 
E 
3 
E 
£ 
31 
5 
E 
8 
E 
c 
2 
E 
o 
[2 


Time since injection (hours) 


Figura 13. Gráfico logarítmico de la actividad de tritio en la bahía (valores porcentuales) en fun- 
ción del tiempo (en horas, con cero correspondiendo al la inyección del trazador). Se ve la gráfica 
de una curva regular que se ajustó a las actividades medidas. 
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2 
Durante la última etapa de desaparición del trazador 7 log, Az (e ) > 0, a diferencia 
Í 


? d? 
de lo que acontece en la primera etapa, durante la cual ri Ár (1)< 0 Recordemos 
Í 


que la actividad total Alt) en el modelo con dos compartimentos clásicos corresponde a la 
actividad A g(£) en el modelo fractal que describe la totalidad del proceso. 


Como podía esperarse, el tiempo de renovación correspondiente al modelo de tres com- 
partimentos resulta mayor que el tiempo de renovación estimado con el modelo conven- 
cional de dos compartimentos, debido al efecto de trampa que el compartimento fractal 
ejerce sobre las moléculas de THO [14]. No obstante, el aumento es significativo pero no 
tan grande como el que se produce al pasar del modelo de un compartimento homogéneo 
al modelo de dos compartimentos homogéneos. Esto se debe a que el efecto de trampa 
comienza a dominar la cinética cuando la mayor parte del trazador ya ha abandonado la 
bahía. 


Efectos de estratificación y de escala sobre los procesos de transporte de ma- 
sa en aguas costeras 


Una vez que sale el agua contaminada de la bahía entre las escolleras Oeste y Sarandí, en 
la época estival encuentra un cuerpo de agua estratificado. Los resultados experimentales 
indican que la mayor parte del agua circula por la capa superior, por encima de la picnó- 
clina. 

Lo contrario ocurre con las aguas residuales descargadas por el emisor de Punta Brava, 
que viajan por debajo de la picnóclina. 

El gradiente vertical en la velocidad horizontal a través de la picnóclina tiende a aumentar 
el transporte turbulento en dirección vertical. Desde el punto de vista cuantitativo este 
aumento se puede resumir estimando una frecuencia de Kelvin-Helmholtz que se calcula 


a partir de las variaciones verticales de las componentes U y V del promedio de la velo- 


or Y (ovY 
cidad horizontal: N=,|| — | +| — 
aa 


La variación brusca de densidad /O a través de la picnóclina tiende a disminuir el trans- 


porte turbulento. Esta disminución puede resumirse estimando una frecuencia de Brint- 


Váisála M= -£(2) (el eje z está dirigido desde el fondo hacia la superficie y 
poz 
£ representa a la aceleración debida a la gravedad). 


A partir de las mediciones en las aguas costeras de la región de Punta Brava se pueden 
hallar valores numéricos aproximados para esas frecuencias en tres regiones: en la capa 
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superior, en la capa inferior y en la picnóclina. Los efectos combinados de ambos proce- 


M?L) 


sos contrapuestos se pueden resumir un número de Richardson local: Ri (z ) = y? ( ) 
z 


Es posible establecer correlaciones (de Mamayev, de Munk y Anderson, etc. [15]) entre 
los promedios, tomados sobre un ciclo de marea, del número de Richardson y del coefi- 


ciente de difusión turbulenta vertical €, (2). A partir de esas correlaciones se puede 


estimar, para los escenarios estudiados, que el valor promedio del coeficiente de difusión 
turbulenta vertical puede disminuir en la interfase (picnóclina) hasta un 80 % de sus valo- 
res en cada una de las capas [5]. En consecuencia la interfase opera como un verdadero 
techo difusivo para el transporte turbulento vertical en la capa inferior (aguas descargadas 
por el emisario) y como un piso difusivo para dicho transporte en la capa superior (aguas 
inyectadas desde la bahía cuando la carrera de marea se dirige hacia el Este). Los resulta- 
dos experimentales relacionados con el transporte de uranina confirman esta conclusión. 


Las primeras etapas de dispersión del trazador, inyectado puntualmente en la capa de 
agua que se encuentra por encima de la picnóclina, se puede describir en una primera 
aproximación, mediante una ecuación aplicable a la concentración clt, Xx, y) de trazador 
promediada sobre la dimensión vertical de esa capa. Si se supone que el proceso de dis- 
persión presenta isotropía horizontal se puede utilizar la ecuación: 


0_ 1.0 g.- 

—clt,r)=-.-—|r-E,lt,r)—clt,r 18 

A) as) 
El coeficiente de dispersión horizontal E, (z, r) que aparece en (19) se suele representar 
así [15]: Ej (t,r)= alt): r? (19) 


El punto correspondiente a r=0se desplaza arrastrado por el movimiento del agua, de 
modo que se trata de un sistema de coordenadas móvil. 

S1 se inyecta una masa M de trazador en r =0 cuando £=0, en una capa de aguas coste- 
ras de espesor h,., la solución de la ecuación del transporte (18) con el coeficiente de 


dispersión dado por (20) es (siendo F(z) la función gama de Euler [16]): 


clt,r)= o . Q-4) oa -expl- o(e).r?-9] Q0) 
, 2 
21 1] —— 
la) 
En (20) se tiene, por definición: o(t)= Q1) 


1 
CY Fale) ar 


Cuando q=0 y alt)= E, es constante, se obtiene un proceso de dispersión gaussiano. 


Cuando q=1 y alt)=P es constante, se obtiene el modelo de Joseph y Sedner que se 


utilizó para la dosificación del trazador fluorescente en los experimentos en aguas coste- 
ras. 
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Las manchas de trazador adoptan formas alargadas e irregulares. Una interpretación posl- 
ble del modelo de mancha con simetría cilíndrica es que describe valores medios sobre un 
conjunto estadístico (un ensemble) de manchas. 

Para aplicar (20) en la práctica se determinan curvas horizontales de iso-concentración, se 
calcula el área encerrada por cada una de ellas, se construyen circulos de igual área y se 
asignan los valores de concentración a los radios de los correspondientes círculos. 


Una vez que la extensión de la mancha es lo suficientemente grande, la dispersión se 
puede describir como un proceso gaussiano, en principio anisótropo. Si v es la velocidad 
horizontal del agua (dirigida paralelamente al eje x y supuesta constante en toda la región 
considerada) y si los coeficientes de dispersión longitudinal (en dirección paralela al mo- 
vimiento advectivo) y transversal son K, y K; respectivamente, entonces: 


EY 2 
al 4-K,:t 4-Kyt 
ho Am te KK 


Esta expresión puede no ser aplicable cuando en la capa de agua considerada hay flujos, 
superpuestos en la vertical, que presentan sentidos opuestos. En ese caso se puede utilizar 
una matriz simétrica de coeficientes de dispersión cuyas direcciones principales ya no 
coinciden con la dirección de la velocidad local y su dirección transversal [17][18]. 


Clt, x, y)= (2) 


Modelo matemático para la estimación del Tso en las aguas costeras próxi- 
mas a Punta Lobos 

Generalmente se supone que la concentración bacteriana c, (1) en un cuerpo de agua 
decae en función del tiempo Í según la fórmula que se deduce del modelo de Phelps [19] 


para los ríos: ¿, (f)= e e *! siendo T,, = = y DÍ ) la dilución. Un trazador con- 
» 
AER (0) 


servativo de concentración c, (e) permitiría medir la dilución: DÍe ) = ] 
colé 


ajustando los parámetros de la fórmula log, Ch (e) =—K,t+108, Cb (0) a los datos expe- 
Cr (e ) Cr (0) 


rimentales, se hallaría el parámetro de auto-depuración K, . Los valores que puede adop- 


Entonces 


tar este parámetro dependen de los procesos que actúen eliminando los microorganismos 
suspendidos en el cuerpo de agua: la absorción de radiación solar y el choque osmótico. 


Debido a las variaciones temporales y espaciales que presentan las aguas costeras de 
Montevideo, la fórmula de Phelps para estimar el Ty, bacteriano no resulta directamente 


aplicable en la práctica. La tasa de dosis de radiación solar que un microorganismo va a 
estar recibiendo (UV y visible) depende de la profundidad a la que se encuentre, de la 
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turbidez del agua, de la estación del año, de la cobertura de nubles y de la hora del día. 
Asimismo, el daño por choque osmótico depende de la distribución de salinidad en el 
cuerpo de agua receptor, de los cambios de posición del micro-organismo arrastrado por 
los movimientos del agua y de su tendencia a sedimentar. 

Por estos motivos se construyó un modelo matemático analítico para describir el transpor- 
te de una población de micro-organismos en un cuerpo de agua estratificado, teniendo en 
cuenta procesos de advección, dispersión, sedimentación e interacciones con el fondo del 
estuario. Se introduce un campo de parámetros locales de autodepuración K, (c, Xx, y, z) 


que depende de las coordenadas cartesianas de posición x,y,z y del tiempo f. Para un 
cuerpo de agua sin piso difusivo (sin una picnóclina marcada), pero con salinidad variable 
en profundidad, luego de introducir algunas simplificaciones (K, (z, Xx, y, z) =K, (z) 
siendo z la coordenada vertical para una columna vertical de agua de profundidad cons- 
tante H), se obtiene la siguiente fórmula asintótica, donde c, (Ea y,z) es la concentra- 


ción bacteriana y C (6%, y, z) es la concentración de trazador conservativo [20]: 
H 


C, (t;x, y,z)dz 
lo l 1% k 2 4 1 Mo, Aso (3) 
g. H > DO * y 08, 57 : A ) 


EXE x, y, 2) dz or Abr 


0 


[K,(2)az 
) a d | 


El coeficiente global de auto-depuración ke. , Tesulta de sumar dos términos. 


a] 24) 


El primero de esos términos es el promedio vertical del parámetro local de auto- 
H 
IE Kp(2): dz 


depuración K, (2) , e 


T 


OIE Th Ñ H es 
En el segundo término —-|1-— | aparece una escala de tiempo 1, =— característica 
lo G W, 


Ss 


del proceso de sedimentación de los microorganismos (w, es una velocidad de sedimen- 


me . . T j 
tación promedio) y el cociente “2 entre el esfuerzo de cizalla 7, sobre el fondo del 


To 


estuario y un valor crítico T, de ese esfuerzo, por encima del cual la erosión con re- 


suspensión se impone al depósito por sedimentación. 
Si hay un piso difusivo la integración se lleva a cabo entre la superficie del cuerpo de 


agua y la cota h de la picnóclina (medida a partir de la superficie). En este caso f, = Ea y 


ws 
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: T : a 
se anula el cociente relacionado con procesos mecánicos que se producen en el fondo 


To 


del estuario. 


Según este modelo, lo que se puede determinar como una estimación de la capacidad de 
auto-depuración a partir de los datos experimentales, es el parámetro global [e que 
comprende dos términos: uno de ellos mide la autodepuración propiamente dicha y el otro 
mide la desaparición o reaparición (según sea 7, menor o mayor que 7.) de microorga- 


nismos en suspensión. Este parámetro permite definir un 7, equivalente así: 


2.3 
Toy = E, 
CONCLUSIONES 


(1) Como consecuencia del brusco aumento en la conductividad (y por tanto de la densi- 
dad) del agua existente en la zona del emisor de Punta Brava cuando se alcanza una de- 
terminada profundidad, durante los meses abarcados por la campaña (verano-otoño tem- 
prano) las aguas residuales vertidas a través del mismo son atrapadas a niveles profundos, 
entre los 3.6 metros (primera detección de la pluma) y los 5.2 metros (quinta detección de 
la pluma). En consecuencia, las playas con una profundidad menor a 3 metros deben estar 
poco expuestas a contaminación proveniente de los caudales descargados a través del 
emisor: la mayor parte de las aguas residuales escapan de la zona de playas a través de 
flujos más profundos. Se comprobó que el emisor sub-acuático operaba de acuerdo a lo 
esperado en las condiciones prevalecientes en el momento del estudio (las usuales en la 
temporada estival) y no sería el principal responsable de la contaminación residual de las 
playas. 


(2) La estimación de los parámetros del nuevo modelo, aplicando un criterio de ajuste 
óptimo a los datos experimentales, presenta una considerable incertidumbre debido a que 
las incertidumbres asociadas a los datos posteriores a las 100 horas a partir de la inyec- 
ción no hace aconsejable utilizarlos en la estimación de parámetros y a que los datos dis- 
ponibles en las primeras 100 horas son relativamente escasos en relación con el número 
de parámetros a estimar. 


(3) De acuerdo con el modelo con un compartimento fractal, la actividad sumada de los 
compartimentos NE y SW, que es lo que se estima a partir de las concentraciones de acti- 
vidad de tritio medidas en los puntos medios de cada celda, disminuye menos que expo- 
nencialmente en función del tiempo. Esta consecuencia del modelo, independiente de la 
incertidumbre en los valores de los parámetros mencionada, suministra una explicación 
posible para las observaciones sobre el comportamiento de la actividad total en la bahía 
durante el período de tiempo comprendido entre 70 y 100 horas de mediciones. 
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(4) La bahía opera como una bomba de descarga intermitente de contaminantes a las 
aguas costeras de Montevideo. De la estimación de los tiempos de transporte de los con- 
taminantes presentes en las aguas costeras provenientes de la bahía (microbiológicos o 
químicos), se desprende que en un lapso no mayor a cuatro días, un agente vertido en las 
dársenas del Puerto de Montevideo alcanza las aguas de la Playa Carrasco. En el caso del 
vibrión colérico, en las condiciones de temperatura y salinidad usuales en la temporada 
estival, cabe esperar una sobre vida promedio de más de cinco días. Debido a casos como 
éste, podría ser necesario tomar medidas sobre el uso de las playas en caso de producirse 
un accidente de contaminación en el puerto. 

Si bien cabe esperar que el incremento en el tiempo promedio de renovación del agua en 
la bahía cuando se pasa del modelo de dos compartimentos clásico al modelo fractal no 
sea muy importante, hay una diferencia importante relacionada con la expresión asintótica 
de la función de densidad de distribución de la probabilidad correspondiente al tiempo 
aleatorio de residencia de una partícula: mientras que en el modelo de compartimentos 
clásico es una función exponencial del tiempo, en el modelo fractal es una exponencial 
dilatada, que tiende a cero más lentamente. Esto puede tener como consecuencia el in- 
cremento en la extensión de la ventana de tiempo asociada a un accidente de contamina- 
ción. 


(5) El proceso de advección-dispersión del trazador, que simula el transporte de un con- 
taminante vertido accidentalmente en forma súbita y localizada en una región con efecto 
de trampa, presenta efectos de histéresis y de escala significativos durante la primera hora 
y media. Luego se puede aplicar un modelo matemático de dispersión gaussiana integrado 
en profundidad. La matriz de coeficientes de dispersión diagonalizada presenta un coefi- 


ciente que toma valores del orden de 1 mi y el otro coeficiente toma valores de un 
s 


orden veinte veces menor. Esta matriz emerge de una estructura de corrientes compleja, 
que puede presentar inversión de flujos en profundidad, cuando se considera una direc- 
ción ortogonal a la línea de la costa, bajo la acción de vientos fuertes y frecuentes prove- 
nientes del S, SE o SW. En este caso los ejes principales de la matriz de dispersión no 
necesariamente coinciden con las direcciones determinadas a partir de la dirección pro- 
medio del movimiento advectivo. Las direcciones de los ejes principales pueden cambiar 
de un punto a otro del cuerpo de agua, si bien de forma continua y en escalas espaciales 
relativamente grandes. Todo esto debe tenerse en cuenta si se va a utilizar un modelo 
analítico de pronóstico, como es el caso de las estimaciones efectuadas con fines de pro- 
tección radiológica. 


(6) El concepto de 7,,, en su acepción original, debe ser refinado considerablemente y 


complementado por otros parámetros para obtener una medida de la capacidad de autode- 
puración aplicable a un cuerpo de agua con una dinámica tan compleja como el estudiado. 
En principio se podrían estimar unos pocos coeficientes cinéticos adicionales, a partir del 
ajuste de los parámetros del modelo matemático presentado en la referencia [20], a los 
datos obtenidos en los trabajos de campo. Para ello se deben utilizar las fases iniciales de 
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los transitorios en la cinética de las concentraciones bacterianas y del trazador conservati- 
vo. 
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El Río de la Plata 


-Posee aproximadamente 320 
km de extensión. 

-Su anchura es variable desde 
varios km en su nacimiento, por 
la confluencia de los ríos Paraná 
y Uruguay, hasta 220 km en su 
desembocadura oceánica. 

-Su profundidad por lo general no 
superan los 20 m. 
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El Río de la Plata: 


Presenta un comportamiento 
hidráulico muy diversificado. 

La influencia de las mareas, los 
caudales fluviales y los vientos 
producen efectos del mismo orden 
de magnitud, pero de importancia 
relativa variable con la zona 
considerada. 


Se pueden distinguir tres zonas: 


Una zona fluvial interna, de agua 
siempre dulce. 

Una zona intermedia con aguas 
predominantemente dulces, con 
estratificación horizontal pero sin 
estratificación vertical. 

Una zona externa de estuario 
propiamente dicho, que presenta 
estratificación tanto vertical como 
horizontal. 


tamento de Montevideo y sus aguas costeras 


-Las aguas costeras de Montevideo se encuentran en la frontera interna de la 
zona exterior, de estuario, donde dicha zona se conecta con la zona intermedia, 
de aguas predominantemente dulces. 


-Las aguas oceánicas, avanzando por el fondo y las aguas continentales, 
avanzando por la superficie, se mezclan influidas por las mareas y los vientos. 


-La franja de aguas costeras de Montevideo, estudiada durante las campañas y 
otros trabajos relatados en la actual presentación, posee unos 5 km de ancho y 
está comprendida entre la boca de la bahía al Oeste y la desembocadura del 
arroyo Carrasco al Este, separados por 18 km de costa. 


-El extremo de la cuña salina formada por agua proveniente del océano se 
desplaza ya sea hacia el Oeste, ya sea hacia el Este respecto de la 
desembocadura de la bahía. 


-Los vientos fuertes del SE empujan la cuña hacia el oeste, pudiendo alcanzar 
la desembocadura del Río Santa Lucía, que forma el límite Oeste del 
Departamento de Montevideo. 


-Cuando los vientos son débiles, la salinidad sigue la propagación de la marea. 


Datos batimétricos de la bahía y de las aguas costeras del 
Departamento de Montevideo próximas a la bahía 
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Tres problemas motivaron los estudios con trazadores realizados en 
las aguas costeras y en la bahía de Montevideo (1) 


La contaminación del agua y los sedimentos en la bahía de Montevideo. 


Cuando se iniciaron los trabajos con trazadores, en 1991, se podía relacionar con: 
-El Puerto de Montevideo. 

-La Refinería de ANCAP y una la Central Térmica de UTE. 

-Dos arroyos que desembocan en la bahía, el Miguelete (21 km de extensión y 113 
km* de cuenca) y el Pantanoso (16 km de extensión y 66 km* de cuenca), cuyas 
aguas y sedimentos se encuentran muy contaminados por desechos originados en 
actividades industriales y en asentamientos de hurgadores a lo largo de sus 
márgenes. 

-Varios colectores de aguas pluviales y cloacales que vertían su contenido a la 
bahía cuando comenzaron los estudios con trazadores, en diciembre de 1991. 


La investigación de la dinámica del agua y los sedimentos en la bahía mediante el 
empleo de trazadores se llevó a cabo con el fin de: 

-Obtener datos utilizables para mejorar la gestión de la calidad ambiental en el 
marco de la implementación del plan de saneamiento de la ciudad. 

-Aportar datos de interés para la autoridad portuaria en relación con la 
implementación del plan de obras civiles para el Puerto de Montevideo. 


Tres problemas motivaron los estudios con trazadores realizados en 
las aguas costeras y en la bahía de Montevideo (2). 


La contaminación de las playas. 


Luego de la puesta en funcionamiento del emisario sub-acuático de Punta Brava 
(Punta Carretas) junto con la eliminación de importantes descargas directas de 
aguas servidas en la costa, a mediados de la década del ochenta, se observó una 
contaminación residual en las playas Ramírez, Pocitos, Malvín y Carrasco. 

La contaminación de esas playas, situadas al Este de la boca de la bahía de 
Montevideo, era significativa y parecía no ser explicable en su totalidad por las 
fuentes de contaminación costera conocidas y cuantificadas. 


Pese a que la concurrencia de los montevideanos a su cinturón de playas decreció 
desde un 78% en 1955 hasta un 54% en 1985, las playas mencionadas 
continuaban y continúan siendo el principal motivo de atracción durante el verano, 
inclusive lo son para muchos turistas. 


Las investigación de la dinámica de la contaminación costera mediante el uso de 
trazadores se llevó a cabo para contribuir a mejorar la implementación de la fase lll 
del plan de saneamiento de la ciudad, incluyendo un segundo emisario sub- 
acuático a instalarse al Oeste de la boca de la bahía, en Punta Lobos. 


Tres problemas motivaron los estudios con trazadores realizados en 
las aguas costeras y en la bahía de Montevideo (3) 


El costo y las consecuencias del dragado del puerto y de los canales. 


El modelado matemático hidrológico y sedimentológico de la bahía y regiones del 
estuario adyacentes, necesario para el plan de obras civiles en el puerto y para la 
eventual reubicación de las zonas de descarga de sedimentos dragados, requiere 
datos de concentraciones de sedimentos en suspensión, tasas de sedimentación y 
movimientos de los sedimentos de fondo, junto con la medición de las velocidades 
del agua, las temperaturas y las conductividades (salinidades) en diferentes 
escenarios. 


El predominio del barro fluido en los sedimentos de fondo conduce a plantearse: 


-La posibilidad de determinar la profundidad náutica sobre la base de la 
identificación de la superficie correspondiente a la concentración límite: 
1.2 tonim 
- La posibilidad de optimizar los ciclos de dragado de las dragas de succión 
extendiendo el tiempo de bombeo, una vez comenzado el reboce, para llevar al 
máximo la masa de sedimentos acumulados en la cántara de la draga, por unidad 
de costo de operación total (incluyendo el recorrido hasta y desde la zona de 
descarga). 


Temas a considerar durante esta presentación (1): 


-Saneamiento y población de Montevideo. Vertidos desde tierra a las aguas de 
la bahía y a las aguas costeras del Río de la Plata próximas a la bahía. 


-Investigación de la dinámica y renovación del agua en la bahía de Montevideo 
utilizando trazadores radioactivos y fluorescentes. 

Identificación de un proceso de transporte de contaminantes en una franja 
estrecha de aguas costeras adyacente a la costa (1991). 


-Investigación de la influencia de los aportes provenientes de la bahía sobre la 
contaminación en las aguas costeras. 

Estudio simultáneo de la dinámica de la pluma de descarga del emisor sub- 
acuático de Punta Brava, también utilizando trazadores radioactivos y 
fluorescentes. 

Identificación de un efecto de trampa o abrigo que presenta la playa Pocitos 
(1992-1993). 


-Estudio del efecto de trampa en la playa Pocitos, empleando trazadores. 
Investigación de las primeras etapas de un episodio de contaminación local, 
simulado mediante el vertido de un trazador. Experimentos para obtener datos 
con el fin de relacionar la micro-hidrodinámica local con la dimensión fractal de 
los bordes de las manchas de trazadores fluorescentes (1994). 


Temas a considerar durante esta presentación (2): 


-Estimación de los parámetros de dispersión de los contaminantes 
conservativos en diversos escenarios representativos de las condiciones más 
frecuentes en las aguas costeras de Montevideo, incluyendo la consideración 
de efectos de escala (1995). 


-Estudios de dinámica de sedimentos en la bahía de Montevideo mediante 
trazadores y fuentes selladas de radiación (1996-1997). 


-Estudios adicionales destinados a obtener información relacionada con 
procesos de dragado, en vista de su posible optimización (1996-1997). 


-Investigación de la capacidad de autodepuración de las aguas costeras para 
un vertido en la ubicación prevista para el emisario sub-acuático de Punta 
Lobos utilizando bacterias y trazadores fluorescentes (2001). 


-Comentarios sobre algunos aspectos de la metodología de protección 
radiológica empleada durante la planificación y la ejecución de las campañas 
de hidrología. 


Saneamiento de Montevideo: Etapas | y Il 


Sistema de saneamiento 


Plt) = 1.8575.000—640.000 xexp| —0.0165977- 6 1975) | (>1975 de tipo mallado unitario 

que en 1993 atendía al 

P(1995)=1.430.000 P(2015)=1.550.000 80% de la población del 
departamento. 


Los desagúes de las 
distintas cuencas del 
sistema se vierten en el 
Río de la Plata. 

Un sistema costero 
recoge los efluentes 
entre el extremo Este de 
la ciudad y Punta 
Carretas, donde una 
planta de bombeo 
impulsa el líquido a 
través de un emisario 
subacuático. 

Un interceptor, en 
construcción durante 
1993, conduce hacia el 
emisario los efluentes 
colectados al Oeste de 
Punta Carretas y hasta la 
ciudad vieja. Los demás 
RIO DE LA PLATA efluentes se vierten en la 
bahía y en los arroyos 
Miguelete, Pantanoso y 
Carrasco. 
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Esquema de vertidos en la zona estudiada 
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Figura 1. Situación de los principales vertidos de aguas residuales 


existentes en la zona costera de Montevideo. 
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Estudios experimentales de la dinámica y la renovación del 
agua en la Bahía de Montevideo (2 al 13 de diciembre de 
1991) 


Mediciones de perfiles de temperatura y conductividad en diferentes estaciones de 
muestreo en la bahía. 


Medida de flujos de agua utilizando trazadores fluorescentes. 
Marcado de la totalidad del agua de la bahía con tritio como trazador. 


Marcado del vertido de aguas servidas del colector denominado12 de diciembre con 
131-1. 


Para el procesamiento de los datos y su representación matemática, la bahía se 
subdividió en celdas a cuyos puntos medios (correspondientes a estaciones de 
muestreo) se les asignaron perfiles de temperatura, conductividad y concentraciones 
de tritio. 


Además, se midieron las concentraciones de tritio en puntos de muestreo distribuidos 
a lo largo del cinturón de playas situadas al este de la boca de la bahía. 


Subdivisión de la bahía en celdas poligonales a partir de las estaciones de control 


situadas en el punto medio de cada celda. A cada celda le corresponde un elemento 
de volumen construido a partir de una profundidad media. 
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Líneas de igual conductividad medida en la capa superficial de la bahía: 
la salinidad como trazador natural para estudiar la mezcla de aguas. 
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Flujos medidos durante una fase de marea creciente. 
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Figura 14. Flujos medidos en la bahía el día 12-XII-91. 


Esbozo del campo de 
Flujos simulados durante una fase de marea creciente velocidades correspondientes a 
una situación sin vientos 
significativos y con una carrera 
de marea creciente. 
El flujo entra al Río de la Plata en 
sentido Oeste con velocidades 
inferiores a 0.5 m/s. 
Se pueden observar el efecto del 
inicio del canal de acceso, 
comprendido entre las escolleras 
Sarandí y Oeste, que modifica la 
dirección del flujo en las regiones 
adyacentes tanto fuera como 
$ dentro de la bahía, así como un 
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Figura 13. Esquema de flujos más frecuente en la bahía durante los 
periodos de marea decreciente. 


Flujos simulados durante una fase de marea decreciente 
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Esbozo del campo de velocidades 
correspondiente a una carrera de 
marea decreciente con vientos 
suaves del SW-S de intensidad de 
aproximadamente 20 km/h. 

El flujo de descarga se dirige hacia 
el exterior del Río de la Plata. 

Es paralelo a la costa en toda la 
zona representada. 

Se pueden observar el efecto, sobre 
el flujo, del tramo Norte-Sur del 
canal de acceso al puerto, debido a 
la salida de una corriente en la boca 
de acceso situada entre las 
escolleras Sarandí y Oeste 
proveniente de la bahía. 

Los resultados se basan en la 
simulación digital de un escenario de 
otoño tardío de 2007 mediante un 
modelo integrado en profundidad 
lación efectuada por el IMFIA 
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Dibujo de la bahía con los principales aportes 
de aguas servidas marcados como puntos 
rojos. 

Se dibujaron las escolleras Sarandí y Oeste. 

Se esbozaron las dársenas del puerto y los canales 
de acceso al puerto y a la refinería de petróleo. 

Se esbozó una posible distribución de direcciones y 
sentidos para el campo de velocidades promedio en 
una capa de 4 m de profundidad. 

La bahía opera como una bomba de descarga 
intermitente de contaminantes a las aguas 
costeras de Montevideo. 
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La dinámica del agua se puede simplificar 


mediante: 

-Un solo compartimento (para estimar un tiempo 
de renovación global) 

-Dos compartimentos, uno noreste (de 
renovación lenta y mezclado directo con vertidos 
de agua dulce de arroyos y colectores de 
saneamiento) y otro suroeste (de renovación 
rápida e interacción directa con el estuario). 

Las características de la batimetría de la bahía y 
adyacencias sugiere la disposición de esos dos 
compartimentos. Las dársenas se dividen entre 
los compartimentos suroeste y noreste. 


Batimetria (m) 
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Variación estimada de la actividad total de tritio presente en la bahía en función del 
tiempo a partir del marcado. La ordenada se presenta en coordenadas logarítmicas. 
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Si se supone aplicable un modelo de un solo compartimento A(t) É 
representativo de la eliminación del tritio, cabe esperar una relación log, — KZ |=-—— 
entre la actividad remanente, la actividad inicial y el tiempo del tipo: A(0) Í, 
A partir de la recta que mejor ajusta a los datos se obtiene una estimación Vv 
del tiempo promedio de permanencia del agua en la bahía. T — 29 horas 
El dato correspondiente a la actividad inicial no se tuvo en cuenta en el es 0 


informe de Viena de 1992. 


Modelos de dos compartimentos para la cinética del tritio en la bahía de Montevideo 


d Ecuaciones de balance de actividad y 


d 
Y, 74, (2) =-0,, a, +0,,:4, di V, (e) = 0, 0,1 +0, volumen para el compartimento NE 


dt 


d d Ecuaciones de balance de actividad y 
v, di A) (+) 7 O, :qa—> O, 24,7 O: 4), dt V, (1) = O, mm O, 2 ES O, e O. 2 volumen para el compartimento SW 


Ecuaciones de balance de 


d 
V()=V,(c)+V,(+) —V()=0, +0,.=0., Á= A, +4, re Y, q; +V, "A actividad y volumen para 


dt : A 
un compartimento único 
V, V, d y 0>+0..) ca Qs 
ESA. AA TA A 
b+ Si p? —4.( Constantes de decaimiento de 
Á, (0) = A(0) Condición inicial para el compartimento NE A, = —————— la concentración de actividad 


en los compartimentos: 


2 
Á, (0) =() Condición inicial para el compartimento SW  b- /b?—4.c La cinética asintótica depende 
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A, de la constante dominante 
a —da ct 
O, A(0) (e7=" —e) 


a, = HH —=:*5—=2 Evolución de la concentración de actividad en el compartimento NE 
YY, (A, M 2,) Simplificación del modelo inicial: 
( A, SRT As ea ) Evolución de la actividad total de trabajamos con los volúmenes promedio de 
Alt) = A(0) itibenTta bahía cada compartimento y suponemos que son 
(A, ll A, ) constantes 
2 Evolución de la actividad total de tritio en la bahía determinada por la constante de 
Alt)= ALO : e %* decaimiento dominante (justificación del procedimiento utilizado en el informe de 
( Ts ) J p 
a Viena de 1992 para estimar un tiempo global de renovación del agua de la bahía) 
La actividad inicial 
d Evolución asintótica de la ” 1 Al) ; aparente que debe 
po Alt) E Alt) actividad total de tritio en la 2 q log, | 4(t) =—— Utilizarse en el modelo 
2 bahía P 1 A, T, asintótico es mayor 
A, quela actividad total 
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Conclusiones principales de la campaña en la bahía (1) 


-La dinámica del agua en la bahía se puede simplificar asumiendo la presencia de 
dos flujos separados, cuyas características dependen de la ubicación variable del 
extremo de la cuña salina en las aguas costeras de Montevideo: 


Un flujo de agua de baja conductividad que ocupa una capa de menos de 3.5 m de 
profundidad. 

Un flujo de agua de conductividad más elevada que ocupa las zonas con 
profundidades mayores a los 4 m. 


Este último flujo penetra a través del canal de acceso al puerto. El otro flujo penetra 
preferentemente a través de la boca Oeste. 


-La estratificación del agua en la bahía parece permanente, lo cual debe ser tomado 
en cuenta al construir modelos matemáticos para la simulación de la dinámica y la 
renovación del agua. 


-Durante la marea decreciente, el flujo que desciende por el estuario, próximo a la 
costa, es desviado por la escollera Oeste y se divide en dos. Se produce entonces 
un flujo secundario que pone en movimiento a una parte significativa del agua de la 
bahía bajo la forma de un flujo circular en el sentido de las agujas del reloj: este 
movimiento es el principal responsable del mezclado y de la renovación del agua. 


Conclusiones principales de la campaña en la bahía (2) 


-En condiciones meteorológicas normales y considerada en conjunto, el agua de 
la bahía se renueva con un tiempo nominal (volumen/flujo) que se halla 
comprendido entre las 15 y las 40 horas, si se asume un único compartimento 
representativo. 


-Si se consideran los compartimentos parciales definidos por las celdas de 
muestreo, el tiempo de renovación disminuye al pasar de las celdas de la parte 
suroeste hacia las celdas de la parte noreste. No obstante, en la mayor parte de 
esas celdas el tiempo de renovación es relativamente corto. 

Las zonas de renovación más lenta corresponden a la desembocadura de los 
arroyos Miguelete y Pantanoso, a las dársenas y a la región comprendida entre 
el espigón F y la central termoeléctrica. 


-El modelo de compartimentos múltiples (uno por cada celda) se puede 
simplificar en un modelo de solo dos compartimentos: suroeste y noreste. 


-Se detectó una señal de tritio en la playa Carrasco (a 18 km de la boca de la 
bahía), claramente por encima del nivel de cuantificación. Este hallazgo 
inesperado sugirió que la bahía podría ser la causante de buena parte de la 
contaminación residual de las playas. 


Aplicación de los resultados al Proyecto de Base Naval del Cerro 


Escollera principal de un nuevo 
puerto (base naval) proyectado en 
el oeste de la boca de la bahía 
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Flujos de volumen 
simulados hacia y 
desde la bahía en el 
estado actual y con la 
escollera proyectada 
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Campañas de marzo y abril de 1993: 


Marcado y estudio de la distribución de las aguas contaminadas que escapan 
de la bahía de Montevideo. 


Propósito principal: contrastar la hipótesis, efectuada luego de analizados los 
resultados de la campaña de diciembre de 1991, según la cual existe un proceso de 
transporte de tipo difusión de rompientes en una franja estrecha a lo largo de la 
costa, capaz de transportar contaminantes desde la desembocadura de la bahía 
hasta las aguas del extremo Este del Departamento de Montevideo. 

Tipo de estudio: euleriano, con cinco transectas fijas, tritio como trazador principal y 
uranina como trazador secundario. 


Marcado y estudio de la pluma de descarga del emisario subacuático de Punta 
Brava (Punta Carretas). 


Propósito principal: estudiar la ubicación y la estabilidad de la pluma de descarga del 
emisario subacuático de Punta Brava, con el fin de valorar su funcionamiento y 
contrastar la hipótesis según la cual dicho emisario no es el causante de la 
contaminación residual observada en las playas de Montevideo, situadas al Este de 
la boca de la bahía, durante la temporada estival. 

Tipo de estudio: lagrangiano, con 131-| como trazador principal. 


Marcación con tritio: dispositivo de rotura de ampollas, toma de 
muestras, pre-concentración, equipo de centelleo para medir 
radiación beta y lanchón de desembarco utilizado. 


' 


Marcado y estudio de la distribución de las aguas contaminadas que escapan 


de la bahía de Montevideo: punto de inyección, transectas de muestreo fijas y 
algunos de los resultados de las mediciones de tritio. 
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Figura I.1. Resultados de los análisis de tritio de muestras medias 


representativas de la totalidad de cada transecta medidas con enriquecimiento 
electrolítico previo. 


Marcado y estudio de la distribución 
las aguas contaminadas que 
escapan de la bahía de Montevideo 
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Ubicación de las cinco transectas, 
con seis puntos de muestreo cada 
una, respecto de las playas y otros 
puntos de interés situados sobre la 


costa. 
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Tabla de resultados de la determinación 
de las concentraciones de tritio en cada 
una de las cinco transectas entre los días 
23 de marzo y 1 de abril de 1993. 


Marcado y estudio de la distribución de 
las aguas contaminadas que escapan de 
la bahía de Montevideo: determinación de 
velocidades en la capa superior mediante 
biplanos. 
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Figura 14. Velocidades del flujo medidas con biplanos 
situados a una profundidad de unos 2,5 metros. 
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Algunos resultados de las mediciones de conductividad y temperatura las 
aguas costeras. 


En las campañas efectuadas en marzo y abril se encontró un cuerpo de agua con una 
estratificación bastante neta. Generalmente se observaron dos capas de agua de diferente 
salinidad, separadas por una picnóclina: una capa superior, de agua más dulce, de río, y una 
capa inferior, de agua más salada, proveniente del océano. 

Al mismo tiempo la conductividad variaba horizontalmente aumentando hacia el Este tanto en 
la superficie como en el fondo y la profundidad de la picnóclina aumentaba hacia el Oeste. 
Una diferencia de conductividad entre el fondo y la superficie de 10 mS/cm a través de una 
columna de agua de 6 m, con una capa de 1 m de espesor a 3 m de profundidad, donde la 
conductividad varía bruscamente en 5 mS/cm (haloclina), es representativa de la situación. 
-Las diferencias de temperatura medidas (menos de 2 * C) entre la superficie y el fondo del 
estuario no inciden en forma significativa en la estratificación, puesto que 2%C corresponden a 
una diferencia de densidad de 0.6 g/l, contra los 7 g/l correspondientes a una diferencia de 
conductividad de 10 mS /cm. 


-De un día a otro se producían variaciones en más o en menos respecto de los valores de 
conductividad obtenidos en el mismo punto y a la misma distancia de la superficie. Estas 
modificaciones, bien conocidas, se deben al avance o al retroceso de masa de agua salada o 
de agua dulce. Pero el patrón de estratificación vertical se mantuvo en todos los casos. 
Durante los días 24 al 28 de mayo se produjo una invasión de agua de baja salinidad 
proveniente del Oeste. La conductividad en la playa Pocitos, durante esos días, no descendió 
a los valores bajos (4 mS/cm) medidos en la playa Ramírez. A partir del día 28 se produjo una 
invasión de agua más salina (conductividades superiores a 35 mS/cm) proveniente del Este. 


Algunos resultados de los experimentos con biplanos y correntómetros 


-Las estimaciones de velocidad efectuadas utilizando biplanos y las mediciones con 
correntómetros muestran en general una diferencia significativa en las velocidades 
promedio entre la capa de agua predominantemente dulce y la capa de agua 
predominantemente salada. Esas diferencias dependen del instante en el ciclo de marea 
astronómica y de los vientos predominantes. En la zona del emisario subacuático , durante 
el segundo experimento con 131-l, las velocidades en la capa superior fueron aumentando 
desde 140 mí/h hasta 310 m/h, mientras que las velocidades en la capa inferior aumentaron 
desde 360 m/h hasta 1200 m/h, para descender luego a 960 m/h hacia el final del 
experimento. 

-El gradiente vertical en la velocidad horizontal a través de la picnóclina tiende a aumentar 
el transporte turbulento en dirección vertical. Desde el punto de vista cuantitativo se lo 
puede resumir estimando una frecuencia de Kelvin-Helmhotz. 

La variación brusca de densidad a través de la picnóclina tiende a disminuir el transporte 
turbulento, lo que puede resumirse estimando una frecuencia de Brunt-Vaisala. 

Los efectos combinados de ambos efectos se pueden resumir en un número a-dimensional, 
el número de Richardson local. 

-Es posible establecer una correlación entre los promedios, tomados sobre un ciclo de 
marea, del número de Richardson y del coeficiente de difusión turbulenta vertical. A partir 
de esa correlación se puede estimar que el valor promedio del coeficiente de difusión 
turbulenta vertical puede disminuir en la interface (picnóclina) hasta un 80 % de sus valores 
en cada una de las capas. En consecuencia la interface opera como un verdadero techo 
difusivo para el transporte turbulento vertical en la capa inferior y como un piso difusivo para 
dicho transporte en la capa superior. Los resultados experimentales relacionados con el 
transporte de uranina confirman esta conclusión. 


Algunos resultados de los ensayos con tritio en aguas costeras de 
Montevideo 


-El agua circula a lo largo de una franja cuyo ancho varía entre 1 y 5 km, con la región más 
estrecha a la altura de Punta Carretas. Los procesos de transporte son muy rápidos a lo 
largo de la costa: en dos o tres días una partícula que parte próxima a la boca de la bahía 
puede alcanzar la playa Carrasco, a más de 15 km de distancia. Un patógeno o un 
contaminante cuya sobrevida promedio en las aguas costeras de Montevideo sea de orden 
superior a la escala de tiempo de esos procesos de transporte, puede alcanzar el límite de 
las aguas costeras de Montevideo con un decaimiento poco significativo. 


-La primera detección de tritio se efectuó a las 2 horas de ocurrida la pleamar siguiente a la 
inyección del trazador. En ese entonces el agua marcada había llegado ya a la transecta 1 
(playa Ramírez) con una concentración elevada en un punto próximo a la costa. Al día 
siguiente no se encuentra tritio en la transecta 2 (playa Pocitos) ni en la zona más próxima 
a la costa de la transecta 3. Pero se halla el trazador en la zona más alejada de la transecta 
3 y en la transecta 4. Esto indica que en la playa Pocitos se produce un efecto de trampa 
que tiene como consecuencia que las descargas locales de contaminantes adquieran más 
importancia allí que en otras zonas de la costa. 


-La aparición de agua de conductividades muy diferentes, marchando en sentidos opuestos 
y mezclándose con el agua marcada no llegan a producir una renovación apreciable del 
agua de la zona costera investigada: el tritio continuó estando distribuido a lo largo de la 
cadena de playas, a concentraciones más bajas pero medibles sin dificultad. 


Marcado y estudio de la pluma de descarga del emisario subacuático de 


Punta Brava: 
Esquema general del emisario 
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Marcado y estudio de la pluma de descarga del emisario subacuático de 
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Colocación del vial con 4.2 Ci de solución de 
i¡oduro de potasio, marcada con 131-1, en el 
recipiente con agua para su posterior inyección en 
la pileta adyacente empleando una bomba 
peristáltica. 


Detección de contaminación externa por el oficial 
de protección radiológica a cargo de supervisar las 
operaciones de marcado del agua del emisario. 


IS 

Detectores de radiación gama utilizados en 
las campañas de hidrología relatadas junto 
con cables y diversos dispositivos de 
procesamiento de las señales. 


Marcado y estudio de la pluma de descarga del emisario subacuático de Punta 
Brava: 
Esquema de los ensayos y resumen de resultados 
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Resumen de una experiencia de marcado de la pluma de descarga del emisario: 


-Se inyectó una solución de ioduro de potasio marcada con 131-| en una de las piletas de 
tratamiento previas al emisario, durante 2 h 45 min, a razón de 0.5 l/min. La concentración 
del trazador era 1377 micro-curio/metro cúbico. 

-La sensibilidad del detector en las condiciones del experimento fue de 2.1 cuentas por 
segundo por metro cúbico y micro-curio. Efectuando un balance de masa se puede estimar 
una concentración de 442 micro-curios por metro cúbico a la salida del emisario, para un flujo 
de 1 metro cúbico por segundo. Esto corresponde a una respuesta de 998 cuentas por 
segundo en el detector. Como el máximo valor medido en el penacho formado a la salida de 
los difusores fue de 181 cuentas por segundo, esto indica que el factor de dilución en ese 
punto era de 5.1. 

-La medición de la actividad del trazador mostró que el agua descargada quedaba atrapada 
por debajo de la picnóclina, en forma completamente estable. Viajaba hacia el Oeste, 
siguiendo la carrera de marea ascendente, aumentando paulatinamente su profundidad y 
expandiéndose horizontalmente más y más. El espesor inicial de la pluma puede estimarse a 
partir del factor de dilución mínimo, del flujo descargado por unidad de longitud del tramo de 
difusores y de la velocidad de la corriente transversal al tramo de los difusores. Si se toma 
una velocidad de 400 m/h (que corresponde al comienzo del experimento), una longitud de 
67 m para el tramo de difusores y el flujo total descargado, se obtiene un espesor de 0.7 m 
para el inicio de la pluma. El espesor estimado a partir de los perfiles de actividad fue de 2 m. 
-Estas dos estimaciones son compatibles por dos motivos: (1) el trazador no viaja en una 
lámina de concentración constante y (2) el detector es capaz de captar fotones gama 


proveniente de puntos relativamente alejados. Aunque el peso del aporte de un elemento de 
volumen al número total de cuentas cae rápidamente con la distancia entre dicho elemento y el detector, 
este último comienza a producir una señal significativamente por encima del ruido de fondo cuando todavía 
se encuentra algo alejado de la nube de trazador. 


Marcado y estudio de la pluma de descarga del emisario subacuático de 
Punta Brava (Punta Carretas): determinación del volumen de respuesta del 
detector de radiación gamma. 
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Figura 3. Variación relativa de la respuesta de un detector 


de centelleo en función del volumen de agua que lo rodea (de 
tector colocado en el centro) para yodo-131. 


El volumen de 
saturación corresponde a unos 500 1. 


Marcado y estudio de la pluma de descarga del emisario subacuático de Punta 
Brava (Punta Carretas): 


Esquema del procedimiento de medición de la actividad en 
función del tiempo y de la posición (tomado de Maggio, Graiño, 


Kopp, Tripoli: “Determinación de los coeficientes de dispersión en el 
Río de la Plata”) 


Trayectoria corregida 


REGISTRO DE ACTIVIDAD 


X 
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Esquema del registro de la actividad del trazador y 
CORRECCION DE TRAYECTORIAS de la construcción de curvas de actividad constante. 


Marcado y estudio de la pluma de descarga del 
emisario subacuático de Punta Brava (Punta Carretas): 
Primera detección de la pluma de descarga. 
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Fig. 8. 
día 31/3/93: primera detección de la pluma. 


Experiencia de marcado con yodo-131 realizada el 


500 m 


Marcado y estudio de la pluma de descarga del emisario 
subacuático de Punta Brava (Punta Carretas): 
Segunda detección de la pluma de descarga. 
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Fig. 9. Experiencia de marcado con yodo-131 0 


realizada el día 31/3/93. Segunda detección 
de la pluma. 


Marcado y estudio de la pluma de descarga del emisario 
subacuático de Punta Brava (Punta Carretas): 
Tercera detección de la pluma de descarga. 
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Fig. 10. Experiencía de marcado con yodo-131 realizada 
el día 31/3/93. Tercera detección de la pluma. 


Marcado y estudio de la pluma de descarga del 
emisario subacuático de Punta Brava (Punta Carretas): 
Cuarta detección de la pluma de descarga. 
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Fig. 11. Experiencia de marcado con yodo-131 realizada 
el día 31/3/93. Cuarta detección de la pluma. 


Marcado y estudio de la pluma de descarga del 


emisario subacuático de Punta Brava (Punta Carretas): 
Quinta detección de la pluma de descarga. 
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0 Fig. 12. Experiencia de marcado con yodo-131 realizada 
el día 31/3/93. Quinta detección de la pluma. 


Conclusiones de los ensayos sobre la pluma de descarga del emisario 
subacuático 


-Como consecuencia del brusco aumento en la conductividad (y por tanto de la densidad) del 
agua existente en la zona del emisario cuando se alcanza una determinada profundidad, 
durante los meses abarcados por la campaña (verano-otoño temprano) las aguas residuales 
vertidas a través del mismo son atrapadas a niveles profundos, entre los 3.6 metros (primera 
detección de la pluma) y los 5.2 metros (quinta detección de la pluma). 

Cuando las velocidades de flujo son bajas, las aguas residuales viajan dentro de una capa de 
menos de 1 metro de espesor que se ensancha progresivamente a medida que el vertido se 
aleja de los difusores. A velocidades mayores aumenta el transporte por difusión turbulenta 
vertical en el fondo del estuario, por lo que se produce un mezclado vertical del agua situada 
entre el fondo y la halóclina (picnóclina). Se detecta experimentalmente a través del perfil 
vertical de actividad del trazador en la 5* detección de la pluma. La dilución del trazador 
aumentó muy rápidamente y las concentraciones alcanzaron en algunos puntos valores por 
debajo de los límites de detección. No obstante, la halóclina no se desvanece: permanece 
como un techo difusivo para las aguas residuales vertidas por el emisario, debido a lo cual la 
transferencia de agua contaminada hacia las capas superiores se produce con dificultad. 

-En consecuencia, las playas con una profundidad menor a 3 metros deben estar poco 
expuestas a contaminación proveniente de los caudales descargados a través del emisario: la 
mayor parte de las aguas residuales escapan de la zona de playas a través de flujos más 
profundos. 

-Se confirma la hipótesis según la cual la descarga del emisario de Punta Brava no sería la 
principal responsable de la contaminación residual de las playas durante la temporada estival. 
-El coeficiente de dilución inicial determinado con trazadores es casi la tercera parte del 
coeficiente estimado mediante la teoría de chorros en medios estratificados. 


Comentarios sobre algunos aspectos de protección radiológica 


-Los trabajos con radioisótopos se llevaron a cabo en respuesta a iniciativas provenientes de 
las autoridades de la Universidad de la República, del gobierno de Montevideo y de la 
autoridad portuaria. 

-En la mayoría de los trabajos con radioisótopos que se describen en esta presentación 
participaron expertos de la sección de hidrología isotópica del OIEA contando con sus propios 
dosímetros. Esto permitió comparar las dosis informadas por el OIEA con las dosis 
determinadas localmente. 

-Para cada trabajo de campo se solicitó y se obtuvo, de parte de la autoridad competente, la 
correspondiente autorización para el uso de radioisótopos sobre la base de la justificación de la 
práctica, comparada con posibles alternativas, su optimización y la limitación de la dosis tanto 
para los trabajadores ocupacionalmente expuestos como para el público en general. 

-Durante las operaciones se aplicaron medidas de seguridad y se realizó un monitoreo 
ambiental, por los oficiales de protección radiológica, de acuerdo con las normas. 

-Se establecieron medidas de seguridad para enfrentar situaciones no planificadas que 
involucraran riesgo de exposición o de contaminación. 

Se produjo un incidente durante las maniobras previas al marcado de la pluma de descarga del 
emisario de Punta Brava con 131-1, que involucró a un funcionario del OIEA. Se resolvió 
satisfactoriamente aplicando las medidas de seguridad previstas para el caso. 

El resto de los trabajos se desarrollaron sin incidentes ni accidentes. 

No obstante, debido a un retardo en el posicionamiento de una de las embarcaciones que 
participaron en el marcado de las aguas costeras con tritio, se efectuó un trasbordo en un lugar 
no planificado. Esto mereció una observación de parte de uno de los oficiales de protección 
radiológica asignados a supervisar los trabajos de campo. 


Campañas de marzo, abril, mayo y junio de 1995: investigación del 
efecto de trampa en la playa Pocitos y determinación de los 
parámetros de dispersión en las aguas costeras del Montevideo. 


Para simular un accidente de contaminación en la playa Pocitos, se inyectaron en 
forma localizada y súbita soluciones de uranina y se midió el transporte advectivo y 
la dispersión de las manchas de colorante. Se repitieron los experimentos en las 
cercanías del emisario subacuático. 

En las campañas llevadas a cabo en la segunda mitad del mese de mayo y en el 
mes de junio, la estratificación vertical típica de la época estival había 
desaparecido, y la diferencia de temperatura entre la superficie y el fondo, de dos o 
tres de grados se había reducido a menos de un grado Celsius. 

Durante la primera etapa del proceso de transporte tanto las características de la 
micro-circulación como la forma inicial de las manchas tenían mucha influencia 
sobre la distribución del trazador (histéresis). El vertimiento de 5 kg de uranina 
disueltos en 200 litros de agua costera durante solamente 5 minutos puede producir 
manchas inicialmente muy alargadas y curvadas, de unos 7 m de ancho y 70 m de 
largo. El efecto de la forma inicial de la mancha puede mantenerse hasta que sus 
dimensiones alcanzan varios cientos de metros. Por este motivo, en relación con 
las posibles consecuencias de un accidente de contaminación localizado en una 
región de aguas costeras con efecto de trampa, presenta cierto interés el estudio de 
la micro-circulación y sus efectos sobre el transporte y distribución del 
contaminante. 


Al principio, se aplicaron técnicas de medición a bordo de las concentraciones de 
trazador, en condiciones de flujo continuo mediante un fluorímetro, durante una 
secuencia de maniobras de navegación efectuadas sobre la mancha en cada una 
de las etapas de evolución seleccionadas. Estas técnicas no permite suficiente 
resolución espacial para estudiar el mencionado efecto de histéresis y detalles 
debidos a los efectos de la micro-circulación. Por este motivo se decidió registrar 
la mancha de trazador, por medios ópticos y desde el aire, para estudiar 
posteriormente la estructura de sus bordes mediante herramientas de análisis de 
patrones. 

Con el propósito de disminuir lo más posible la histéresis, se diseñó y construyó 
un corral para la inyección de trazadores, que resultó muy eficaz, por lo cual fue 
utilizado en las campañas siguientes. 
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Corral para inyectar trazadores de modo de minimizar efectos de histéresis y 
facilitar la investigación de los efectos de escala en aguas costeras. 


Para determinar los parámetros de la dispersión se efectuaron dos inyecciones en 
la capa superior del cuerpo de agua estratificado y una inyección en el cuerpo no 
estratificado. Se encontraron efectos de escala análogos a los encontrados en 
otras zonas del Río de la Plata y en otras aguas costeras. 
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Representación de líneas de igual 
concentración de trazador 
fluorescente determinada a bordo 
por medidas en flujo continuo. 


Se describió el efecto de escala suponiendo un 
coeficiente de dispersión horizontal del tipo E, = 4:r? 
donde q representa el parámetro de efecto de escala. 
Aplicando el método de los círculos equivalentes, se 
obtuvieron valores cercanos a 1 al comienzo del 
proceso. Los valores fueron disminuyendo 
continuamente hasta anularse cuando la nube de 
trazador se había hecho lo bastante extensa, 
aproximadamente luego de transcurrida una hora y 
media desde el momento de la inyección. A partir de 
esa etapa se podía asumir un proceso de gradiente. 
Para los coeficientes de dispersión longitudinal se 
obtuvieron valores próximos a 1 mó) K, 0.83 m/ ) y 
cerca de treinta veces menores AN coeficientes 
de dispersión transversal (X7 =0.018%/. ), 
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Algunas conclusiones sobre las campañas para determinar los 
parámetros de dispersión en las aguas costeras de Montevideo: 


-El proceso de advección-dispersión del trazador, que simula el transporte de 
un contaminante vertido accidentalmente en forma súbita y localizada en una 
región con efecto de trampa, presenta efectos de histéresis y de escala 
significativos durante la primera hora y media. Luego se puede aplicar un 
modelo matemático de dispersión gaussiana integrado en profundidad. 


-La matriz de coeficientes de dispersión diagonalizada presenta un coeficiente 
que toma valores del orden de ] m/ y el otro coeficiente toma valores de un 
orden veinte veces menor. ñ 


-Esta matriz emerge de una estructura de corrientes compleja, que puede 
presentar inversión de flujos en profundidad cuando se considera una 
dirección ortogonal a la línea de la costa bajo la acción de vientos fuertes y 
frecuentes provenientes del S, SE o SW. En este caso los ejes principales de 
la matriz de dispersión no necesariamente coinciden con las direcciones 
determinadas a partir de la dirección promedio del movimiento advectivo. Las 
direcciones de los ejes principales pueden cambiar de un punto a otro del 
cuerpo de agua, si bien de forma continua y en escalas espaciales 
relativamente grandes. Todo esto debe tenerse en cuenta si se va a utilizar un 
modelo analítico de pronóstico, como es el caso de las estimaciones 
efectuadas con fines de protección radiológica. 


Campañas de estudio de sedimentos 
en la bahía de Montevideo y las zonas 
de descarga del material dragado. 


(1) Seguimiento de material dragado 
marcado y descargado por la draga. 


Determinación de parámetros de dispersión de la 
nube de partículas en suspensión, estimación de 
tasas de sedimentación y estudio del movimiento 
del material sedimentado. 


(2) Determinación de la profundidad náutica mediante fuentes selladas de 
radiación en situaciones en las que el sonar no resulta aplicable. 
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La capa de barro fluido en el fondo de los canales de navegación en la Bahía de 
Montevideo posee características que justifican el empleo de una sonda de 
densidad para la determinación de la profundidad náutica con el fin de ajustar la 
profundidad de dragado teniendo en cuenta el calado previsto de los barcos. 


(3) Optimización del tiempo de bombeo en dragas de succión. 


Composición promedio del 
material que ocupa las zonas 
a dragar en la bahía de 
Montevideo: 


Arena: 5% 

Limo: 80% 

Arcilla: 15% 

Arena + limo grueso: 25% 


Se justifica determinar la 
curva de acumulación de 
material en la cántara de las 
dragas que operan por 
succión continua con el fin 
de estimar el tiempo óptimo 
de operación de las bombas 
y minimizar el costo de cada 
ciclo de dragado. 


Supongamos que el costo de una operación de dragado es proporcional a la duración de los 
procesos que la componen. Consideremos las rectas que cortan a la curva de acumulación de 
masa en función del tiempo y cortan el eje del tiempo en un instante negativo cuyo valor 
absoluto representa el tiempo de operación. Entonces la recta tangente define el tiempo 
óptimo de operación de las bombas, porque se obtiene un valor máximo para el cociente entre 
la masa de sedimentos acumulados y el costo del dragado. 


Investigación de la capacidad de autodepuración de las aguas 
costeras de Montevideo en la región de Punta Lobos mediante la 
determinación delZop (junio y julio de 2001). 


El objetivo de los trabajos fue la estimación del tiempo de inactivación ( 194 ) de una población 
bacteriana (E. coli): 
To = 2 ¿(o cel0) ¿ku Dolo ero) og. 2 a, 24109, 20) 

K DÍt) Er (+) cy (1) cy (0) 
Se llevaron a cabo cuatro campañas durante los mese de junio y julio de 2001. En todos los 
casos se inyectó, en la zona prevista para el emisario subacuático, en forma súbita y puntual, 
una masa de efluente conteniendo la población bacteriana, acompañada de un trazador 
fluorescente y conservativo (rodamina WT) para determinar la dilución por el proceso de 
advección-dispersión. Se utilizó uranina como trazador secundario cualitativo, para visualizar la 
ubicación del vertimiento efectuando marcados sucesivos al finalizar cada etapa de 
mediciones. En cada etapa se tomaron muestras para determinar unidades formadoras de 
colonias y concentración del trazador en laboratorio. Además se determinó la posición de cada 
punto de toma de muestras, los perfiles verticales de velocidad de corriente, de temperatura y 
de conductividad, la velocidad local del viento, la turbidez del cuerpo de agua, el pH, la DBO y 
otras variables características de cada escenario en el que se produce la advección, la 
dispersión y la inactivación bacteriana. 
Según la profundidad a la que se encuentre un micro-organismo, la turbidez del agua, la 
estación del año, la cobertura de nubles y la hora del día, va a ser la tasa de dosis de radiación 
solar que va a estar recibiendo (UV y visible). Asimismo, el daño por choque osmótico depende 
de la distribución de salinidad en el cuerpo de agua receptor, de los cambios de posición del 
micro-organismo arrastrado por los movimientos del agua y de su tendencia a sedimentar. 


Debido a las variaciones que presentan las aguas costeras de Montevideo, la fórmula para 
estimar el Lo bacteriano no resulta aplicable en la práctica. 

Por este motivo se construyó un modelo matemático para describir el transporte de una 
población de micro-organismos en un cuerpo de agua estratificado, teniendo en cuenta procesos 
de advección, dispersión, sedimentación e interacciones con el fondo del estuario. 

El modelo introduce un parámetro local de autodepuración K, o que depende delas 
coordenadas de posición y del tiempo. Para un cuerpo de agua sin piso difusivo debido a una 
halóclina muy marcada, pero con salinidad variable en profundidad, se obtiene, luego de 
introdu ir algunas simplificaciones, la siguiente fórmula asintótica: 


EXé ld ñ 
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Si hay un piso difusivo la integración se lleva a cabo entre la superficie y la cota de la halóclina. 
Además, en el término relacionado con la sedimentación no aparecen el cociente entre el 
esfuerzo de cizalla sobre el fondo del estuario y el valor crítico de ese esfuerzo, por encima del 
cual la erosión con re-suspensión se impone al depósito por sedimentación. Con fines de 
orientación, para el diseño de los experimentos de campo y para la estimación del orden de 
magnitud del tiempo de inactivación en diferentes escenarios durante el transcurso de un año se 
utilizó la expresión: am A S 


a, e (1-(0.65,0?). 1-2 
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_Datos Fluorimétricos obtenidos en el Laboratorio de la DINATEN a partir de 
muestras tomadas por la 1.M.M. 


Cuarta campaña en la región de emplazamiento prevista para 


el emisario de Punta Lobos: mediciones de bacterias y trazador cónicos: 1- supericiar (20 cn pa , 
2-1.5m 0.303ppb 
BACTERIAS 


C. FECALES A dilución 1:200 
ufc/100mL - AA 34 55' 23" 5 ; dilución 1:200 
561605" w [—a000 | diución1:200 
Superficie 4451 


S34":55"23” 1,5m ¿980 | dilución 1:100 
W 56” 16' 09" 3,5m 34 55' 20" S y dilución 1:100 
columna E 3: 56 16'05" W 
Superficie 
S 340 55' 20" 1,5m —. $ dilución 1:100 | 3.957 
cd ri 56 16' 02" W 
columna 
Superficie 
S 349 55' 19" 1,5m 
W 56" 16' 02" 3,5m 
columna 
Superficie | j dilución 1:100 
S 34? 35' 18" 1,5m ¿ l 34 55' 18" S 
W 56' 15' 58" 3,5m | 56 15' 51" W 
columna 
Superficie | 
S 340 55' 18" 15m | 34 55 18" S 
W 56 15 51" 3,5m 56 15 46" W 
columna 
Superficie | Ea 
S 340 55'18" 15m | be e da Y 
W 56" 15' 46" 3,5m l 
columna 
Superficie | notorio color de 
S:3A7D5' 17" 1,5m 34 55' 18" S : ¡uranina 
W 56" 15' 40" 3,5m 56 15' 30" W 
columna 


34 55' 18" S 
56 15' 58" W 


Superficie 
S 34" 55' 18" 1,5m 5 34 55' 18" S 
W 56* 15' 30" 3,5m 56 15' 18" W 
columna | 


Superficie 

S 34" 55' 18" 1,5m 

W 56 15' 18" 3,5m 
columna 
Superficie 

5345516" 15m 

W 56" 14' 51" 3,5m 
columna 
Superficie 

S 34" 55' 16" 13m 

W 56" 14' 23" 3,5m 
columna 


34 55' 16" S 
56 14' 51" W 


34 55' 16" S 
56 14' 23" W 


34 55 14" S 
56 13' 58" W 


Superficie 


Yi id 2 TANQUE [nuestra DNTN_ [286 lución 1:100.000 
W 56" 13' 58 3,5m (0,1mi en1000, y 
columna 200L ¡muestra IMM 31,9  !10mlen 100) 


Medidas de Temperatura, Conductividad y Salinidad: 


PROF. (m) 


superficie 


SALINIDAD 


TEMP.(*C) 


CONDUC. (mS/cm) 


superficie 


0,5 


COORDENADAS 


Lat./Lon.: 
34%55' 23" S 
56" 16' 05" W 
Anemómetro: 1,38 m/s. 
Hora 10:40 


Lat./Lon.: 
34055' 199" S 
56% 15' 273" W 
Anemómetro: 2,21 m/s. 
Hora 10:40 


Lat./Lon.: 
34955' 179" S 
56" 15' 097" W 

Hora 13:25 


Lat./Lon.: 
34055" 163" S 
56" 14' 580" W 
Anemómetro: 2,47 m/s. 
Hora 13:35 


Medidas de Corriente 


PROF. (m) 


superficie 


TEMP. (*C) 


VELOCIDAD (cms) 


DIRECCIÓN 


11,25 


11,1 


11,05 


11 


2,2 


10,8 
84 
9/2 


274,3 

118,1 

126,7 
188 


— == 00 A 
10,8 249,9 | 


260,3 
274,3 
274,3 


COORDENADAS 


34* 55' 23" S 
56" 16' 05" W 
Hora 10:45 


34" 55' 229" S 
56" 16' 09" W 
Hora 11:30 


34% 55' 199" S 
56" 15' 273" W 
Hora 12:45 


Algunas conclusiones sobre la determinación de la capacidad de 
autodepuración de las aguas costeras en la región de Punta Lobos. 


-Se estimaron valores para el Ese comprendidos entre 6 y 12 horas, tomando 
valores mínimos en los escenarios correspondientes a las dos últimas campañas 
y valores máximos en los escenarios correspondientes a las dos primeras. 


-Si bien las temperaturas del aire durante las últimas dos campañas eran 
relativamente elevadas para la época del año, la temperatura del agua no lo 
reflejaba, como cabía esperar. 


-En las primeras dos campañas la turbidez del cuerpo de agua era 
significativamente mayor que la turbidez durante los dos últimos trabajos de 
campo. 


-Con una metodología análoga a la utilizada en análisis compartimental, se 
podrían estimar unos pocos coeficientes cinéticos adicionales a partir del ajuste 
de los parámetros del modelo matemático a los datos obtenidos en los trabajos 
de campo. Para ello se deben utilizar la fase inicial de los transitorios en la 
cinética de las concentraciones bacterianas y del trazador conservativo. 


Fin de la presentación 


¡Muchas gracias! 
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1-RESUMEN 


El trabajo presenta los resultados obtenidos, a través de una 
experiencia realizada mediante la inyección de trazadores (uranina 
y tritio) en la desembocadura Este de la Bahia de Montevideo. 


ss El objetivo fue definir si existe O no influencia de la 
contaminación de las aguas de la bahía sobre las playas de la zona 
balnearia de la ciudad localizada al Este. 


Los resultados obtenidos indican que efectivamente los 
contaminantes que se encuentran en la bahía, se distribuyen a lo 
largo de la costa en un fenómeno de dispersion paralelo a la costa, 
a poca distancia de la línea de rompiente, para las condiciones 
climáticas vigentes, que fueron las promediales de la temporada 
estival. 


Se describen los fenómenos de contaminación que motivaron el 
estudio, los fenómenos de transporte más relevantes de la costa de 
Montevideo, la metodología experimental utilizada (de muestreo y 
laboratorio) y el soporte matemático de las hipótesis manejadas. 


NOTA: La presente investigación experimental se llevó a cabo, 
a pedido de la Facultad de Ingenieria de la Universidad 
de la República en el marco del Proyecto uURUTg1007 que 
la Dirección Nacional de Tecnología Nuclear lleva a 


cabo con la asistencia de la Agencia Internacional de 
Energía Atómica. 


2-OBIJETO 


El objetivo de este estudio p es determinar si la contaminación 
de las aguas de la Bahía de Montevideo, afectan las Playas 
localizadas al Este de su desembocadura, mediante la utilización de 
trazádores adecuados. 


3-DESCRIPCION DE LA PROBLEMATICA 


La ciudad de Montevideo, capital de la República Oriental del 
Uruguay está situada al Sur del país con una localización 
geográfica aproximada de 35 grados latitud Sur y 56 grados longitud 
Oeste . 

La ciudad limita al Sur con el Estuario del Río de la Plata, la 
costa posee una bahia con puerto natural en su extremo Sur Este y 
una cadena de hermosas playas que representan uno de sus recursos 
naturales más importantes, desde el punto de vista de la industria 
turística y como ¿ona de recreación y esparcimiento por excelencia 
de la población durante la temporada estival. (Ver fig.1) 

La ciudad posee un sistema de saneamiento de tipo unitario que 
atiende en la actualidad a 1:000.000 de personas lo que 
representa aproximadamente el 80% de la población del 
departamento. 


Los desagúes de las distintas cuencas Se vierten al Río de la Plata 
de diversos formas. la primera de ellas y la más importante es el 
sistema costero, recientemente inaugurado, que recoge los efluentes 
entre el extremo Este de la ciudad y Punta Carretas donde una 
planta de mbeo impulsa el líquido a través de un emisario 
«subacuático de 2.300 m de longitud a su punto de disposición final. 
Al Oeste de Punta Carretar , un ¡interceptor actualmente en 
construcción deberá conducir al mismo punto de disposición los 
efluentes colectados desde el borde de la bahía de Montevideo 
(Calle Guarani,en el comienzo de la Escollera de cierre del 
Puerto). 

Debe destacarse que los trabajefewes de saneamiento de la costa, 
asi como los diseños realizados, han tenido como finalidad la 
descontaminación bacteriana de las aguas de recreo, razón por la 
cual no se enfatizaron los estudios en relación a la incidencia de 
contaminantes conservativos en las aguas costeras. 


El resto de los líquidos generados en la ciudad se vierten,un 
pequeño porcentaje al arroyo Carrasco en el extremo Este de la 
costa, 
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y el saldo tienexs como destino final la bahía de Montevideo, 


Festa es un espejo de agua semicerrada de 1200 Ha de extensión, Cuya 
profundidad oscila entre 0,50 y 3.50 ma excepción de los canales 
del Puerto de Montevideo, donde mediante dragado y atendiendo a 
las necesidades de la navegación de ultramar se llega a 10 m 
aproximadamente. 


El desarrollo urbano de la ciudad se dio históricamente alrededor 
de este cuerpo de agua, cuya principal característica es la de 
ser la estructura geográfica del Puerto, centro logístico del 
desarrollo mercantil del país. 

Las instalaciones portuariasf muelles, depósitos, playas de 
contenedores, vias de comunicación etc.) abarcan un barrio: La 
Aguada. Actualmente Autoridades portuarias están remodelando 
dichas instalaciones con vistas a agilitar ééer la tarea 
portuaria en su conjunto. 


Naturalmente la intensa actividad del puerto promovió en las 
proximidades del mismo el surgimiento de servicios relacionados, 
como por ejenplo astilleros, dique seco, barracas, etc. 


A orillas de la bahía se encuentra además e la refinería de 
ANCAP, que posee un muelle de atraco para pequeños buques tanque 
destinados al abastecimiento de derivados de petróleo al litoral 
del país. Dicho muelle se comunica con el Puerto a través de un 
canal de navegación que atraviesa la bahía. 


Los usos actuales del cuerpo de agua, además de los portuarios 
son los siguientes: - suministro de agua para enfriamiento (la 
Central Batlle de UTE y la refineria de ANCAP utilizan el agua de 
la bahía para esta finalidad), pesca artesanal y deportiva, 
navegación deportiva. : : 

La única zona de baños actualmente en uso es la Playa del 
Cerro, situada en el extremo Oeste, pues el resto del espejo de 
agua no posee caracteristicas de balneabilidad, debido no solamente 
a su alto nivel de contenido bacteriano sino además por la 
presencia de cantidades muy importantes de residuos sólidos, 
aceites, grasa e hidrocarburos. 


Las infraestructuras de abastecimiento de agua y saneamiento 
acompañaron el crecimiento de la cuidad con carencias históricas 
sobre todo en Jo relacionado con el segundo rubro, situación que 
al agravarse llevó a la bahía al actual estado de deterioro 
ambiental generalizado. 

En partícular, el saneamiento de una parte muy ¡importante de 
log barrios circundantes, Ciudad Vieja (incluido el Puerto), 
Centro, Arroyo Seco, fue realizado a fines del siglo pasado y 
principios de este. 

La red construida, de tipo mallado unitario, vertía los 
efluentes sin mayores cuidados a la bahía y a la costa Sur de la 
ciudad. 

Como parte de las obras de implantación del Puerto, se 
construyeron a principios de siglo dos colectores con la finalidad 
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áe eliminar los vertimientos de líquidos cloacales a la zona 
portuaria pues se los considaraba responsables por algunos de los 
problemas de salud pública que por ese entonces asolaron la ciudad 


El primero de los citados fue construido en forma perimetral 
a la península de la cludad vieja, llevando los afluentes a la 
altura de la confluencia de la costa Sur con la calle Guaraní. El 
segundo partiendo del centro de la ciudad en la Avenida Rondeau, 
atraviesa la parte alta de la cuchilla y desagua en la rambla Sur 
a la altura de la calle Paraguay. 

El sistema construido cumplió con el cometido de alejar los 
efluentes de la bahía hasta mediados del siglo; a partir de esa 
época la situación se fue deteriorando en forma paulatina y 
sostenida. 


Este escenario se fue consolidando debido a Jos aumentos de 
caudal provocado por la extensión del problema y se agrava POr el 
fenómeno de la sobrecarga de los colectores debido el aumento de 
la población atendida. 

Cuando se llegó al límite de capacidad del sistema se 
pretendió entonces subsanar el problema hidráulico con vertimientos 
de alivio a la bahía, aumentando de esta forma la cantidad de 
contaminantes vertidos a sus aguas. (Ver fig. 2) 


A esta situación relacionada con la red de saneamiento de la 
ciudad (correspondientes a las denominadas cuencas Paraguay alto y 
Ciudad Vieja), se le suman los vertimientos del barrio del Cerro 
(de escasa importancia) y de los dos principales cuerpos de agua 
interiores del Departamento, los arroyos Miguelete y Pantanoso. 


El primero aporta principalmente vertimientos de aguas 
cloacales de origen doméstico y una cantidad de residuos sólidos 
muy importantes povenientes de los barrios de  hurgadores 
establecidos a sus orillas en el punto medío de su cauce. 


El segundo aporta una capacidad muy grande de efluentes de 
origen doméstico AÑ e industrial, debido a la concentración de 
industrias instaladas sobre todo en su margen izquierdo, ye destaca 
que además de la materia orgánicas se encuentran presentes otros 
contaminantes de mayor riesgo como por ejemplo Cromo aportado: por 
las numerosas curtiembres asentadas en su cuenca. En resumen se 
puede decir que la cantidad de materia orgánica (medida en DBO 
5,20) vertida en la actualidad a la bahía se estima en 39 Ton/día 
por el sistema de vertedores de alivio de la cuenca Paraguay Alto 
y Ciudad Vieja, 19 Ton/día Arroyo Miguelete, 15 Ton/día Arroyo 
Pantanoso y 4 Ton/día barrio Centro. 

En relación a contaminantes químicos como por ejenplo el Cromo 
se estima que se vierten a la bahía 1.000 Kg/día, de los cuales el 
80% como jones solubles. 


El cuerpo de agua comienza a dar signos de sobrecarga en 
relación a su poder de recuperación natural, apareciendo problemas 
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á 4 
de olores por disminución de oxígeno, sobre todo despues de las 
iluvias cuando el aporte de los arroyos se hace mayor y se acentúa 
la sedimentación de materia orgánica. 


A modo de conclusión puede decirse que el nivel de 
contaminación de la bahía, debido a la sobrecarga de materia 
orgánica y de elementos químicos que se vierten en ella, provocan 
la degradación del cuerpo de agua con la aparición de olores 
ofensivos en algunos puntos. Esto sumado a la pésima condición de 
la costa provotada por la presencia de residuos sólidos, aceites 
y grasas lo convierte en un foco de contaminación para el resto de 
la costa del Departamento. 


En particular, durante un estudio de la dinámica del agua en 
la bahía de Montevideo, efectuado a fines de 1991, se encontró una 
señal de tritio (trazador empleado para marcar el agua de la bahía 
) en el agua de la orilla de la playa Carrasco, á 15 Km de la boca 
de la bahía, antes de los tres días posteriores al marcado. (Plata, 
Kurucz, Suárez Antola y otros, 1992). Este hecho, así Como 
resultados de otros estudios preliminares (incluyendo simulación 
digital del transporte) sugerían la existencia de una difusión de 
rompientes en una franja estrecha a lo largo de la costa, capaz de 
transportar contaminantes desde la desembocadura de la bahía en el 
Rio de la Plata hasta las aguas del extremo este del Departamento 
de Montevideo. 


El presente trabajo tuvo Como finalidad confirmar esta 
hipótesis para lo cual se desarrollaron las tareas 
interdisciplinarias e interinstitucionales de campo, de laboratorio 
y de gabinetes que se describen a continuación. 


4-METODOLOGIÍA 


Para la realización de los trabajos de campo correspondientes, 
se utilizaron dos tipos de trazadores, uno químico: uranina y otro 
radiactivo: tritio. Las ventajas del tritio como trazador para 
este tipo de estudios son su no degradabilidad en períodos 
relativamente largos, compatibles con el tiempo de la experiencia 
y como aparece formando parte de una molécula de agua, el 
comportamiento de las moléculas de agua marcadás és idéntico al de 
las moléculas no marcadas. En cuanto a la uranina tiene la 
desventaja de degradarse en el tiempo bajo la acción de varios 
factores, principalmente la acción de la luz. Si se tratara de 
marcar un gran volumen de agua, caso de un estuario, la cantidad 
necesaria de esta sustancia seria muy grande, no justificandose su 
utilización debido a su costo. Teniendo en cuenta la sensibilidad 
del equipo de medida, se calculó necesario 12 Kg de uranina para la 
experiencia. 

Interesa poder establecer, en caso de existir, un vínculo 
entre los comportamientos de ambos trazadores. En efecto, el 
seguimiento de los procesos de advección y dispersión de un 
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trazador directamente integrado a la molécula de agua (tritio) 
suministra información directa para evaluar el transporte de 
elementos conservativos y un marco de partida para evaluar el 
transporte y la degradación de elementos no conservativos. La 
uranina en conjunto con el tritio, permite una comparación cruzada 
de interés para planificar futuros estudios experimentales 
complementarios. 

Fueron utilizados 100 curios de agua tritiada y 12Kg uranina 
para el marcado 

El procedimiento utilizado para el vertimiento de la uranina 
en el cuerpo receptor, fue el tradicional de diluirla previamente 
y luego verterla, operación que no presenta mayores dificultades 
operativas. En relación al tritio hubo que diseñar un dispositivo 
especial que denominamos "Rompedor de ampollas" (Ver fig.3), 
formádo por dos caños de hierro (uno dentro del otro) de 3m de 
largo, finalizando en un depósito donde se coloca la ampolla a 
romper. Este sistema fue diseñado como un dispositivo de acción 
segura. ' 


El punto de inyección de ambos trazadores, se seleccionó 
teniendo en cuenta los resultados obtenidos en la primera campaña 
del proyecto "Isotopos en Hidrología".URU/8/007 del OIEA. En dicha 
campaña se observó una alta renovación de las aguas de la banía (48 
horas). El punto de inyección esta ubicado a unos 500m de la boca 
de la bahía y a 100m al este del canal de acceso al puerto (se 
quiso evitar de esta forma la influencia del canal que provocaría 
perdidas de trazador). 


Se determinaron cinco transectas de 3 km de longitud en las 
playas: Ramirez, Pocitos, Honda, Verde y Carrasco; en cada una de 
ellas se distribuyeron en forma proporcional seis puntos de 
muestreo. (Ver fig. 4) 


En cada transecta se colocó una boya en el extremo más alejado 
de la costa, desde ella y en forma perpendicular a la costa a 
distancias predeterminadas, se encontraban los demás puntos de 
muestreo. 


La localización de la boya y demás puntos se realizaron 
mediante dos unidades de GPS con un error máximo de 30m 


La metodología empleada para la extracción de muestras, debido 
que se trataba de 30 puntos y de una distancia total de más de 15 
km, fue la de integrar en cada punto en profundidad. 

Para ésto se decidió utilizar una manguera plástica lastrada 
en su extremo inferior, cuyo diámetro interno se calculó a través 
del siguiente procedimiento. 

La tensión superficial generada en tubos de pequeños diámetros 
se relaciona con la altura del ascenso capilar h a través de la 
formula: 
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e ES ” 
Donde: 6, tensión superficial, para agua limpia y temperatura 
normal =0,071 N/m. 


a, ángulo de contacto que representa la adhesión entre 
el agua y el tubo, a=20*. 


d, diámetro del tubo, a determinar. 


h, altura capilar. (Se consideró la mayor profundidad 
de muestreo h=5m). 


t, peso unitario del agua. 


para estos valores se obtuvo un diámetro necesario, de5,3m.m. 

El diámetro mínimo encontrado en el mercado fue de 10m.m, para 
el cual se obtiene, hmax =2,6m. Para Jos puntos de mayor 
profundidad, las perdidas de agua eran excesivas. Fue necesario 
entonces modifica el método de extracción de muestras: se obstruyó 
parcialmente el extremo inferior del tubo, reduciendo su 
diámetro, evitando fugas importantes de agua, y se colocó una cuerda 
adicional atáda del extremo inferior del tubo plástico, con la 
finalidad de una vez sumergido este y tapado su. extremo superior, 
tirar de ella levantado así la columna en su totalidad. 


5-TRABAJO DE CAMPO 


Tratándose de más de 15km de costa a recorrer desde la primer 
transecta a la última, el total de muestras obtenibles por día 
dependía fundamentalmente de la embarcación disponible. Se contó 
básicamente con dos embarcaciones de características muy 
diferentes, una lancha de desembarco muy cómodo para trabajar, pero 
lenta, ( se tardaba aproximadamente de 8 a 9 horas en extraer todas 
las muestras). Esta embarcación fue muy útil para el vertido de 
ambos trazadores por su seguridad y estabilidad y para los primeros 
días de muestreo cuando las distancias a recorrer eran cortas. La 
otra embarcación, era un pequeño yate de condiciones óptimas para 
este tipo de trabajo, pues desarrolla velocidades altas, lo cual 
permitió reducir casi a la mitad el tiempo de muestreo. 


El marcado con uranina se realizó diluyéndola en un tanque de 
200 litros y utilizando una manguera para su distribución. Cabe 
destacar que como se utilizó la uranina en forma de polvo y esta es 
muy volátil, se corrió peligro de contaminar los equipos, 
frascos,etc, pudiendo interferir esta situación falseando las 
modidas. 


La ampolla de 100Ci de agua tritiada se colocó en el rompedor, 
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este se sumengió aproximadamente 2m y se procedió a romperla. 


El primer día se muestreó solo la primer transecta situada a 
5 Km del punto de inyección. El segundo día se muestreo hasta la 
segunda transecta. En ambos días se encontró uranina, se decidió 
entonces muestrear todas lastransectas a partir de tercer día, 
debido a su rápido desplazamiento paralelo a la costa. 


En cada punto de muestreó se realizaron las siguientes 
tareas: 1) extracción de muestras integrada en profundidad 
2) lectura de posición, GPS 
3)medidas de conductividad y temperatura 
1 volumen extraído en 1), se vertía en un recipiente, se 
mezclaba y se llenaban dos frascos de 50m1 (uno para medición de 
uranína, otro para medición de tritio). En cada frasco se colocó un 
código donde se registraba hora, número de transecta, punto dentro 
de cada transecta y si era integrada en profundidad o puntual. 


En algunas zonas donde se observd el mezclado de aguas de 
distinta conductividad, se extrajeron muestras de superficie de 
ambos lados y en la mezcla misma, para identificar el 
comportamiento de los trazadores. 


6-TRABAJO DE LABORATORIO 


Las muestras de uranina debidamente protegidas de la 1uz solar 
desde su extracción, fueron llevadas al laboratorio para su 
posterior análisis en el Fluorímetro. La ventaja de este trazador 
frente al tritio, es su rápida determinación. En efecto, algunas 
horas después de llegada la muestra al laboratorio se pueden tener 
resultados confiables de las medidas. Además, aunque en condiciones 
más defavorables desde el punto de vista analítico, se puede 
realizar mediciones in situ a tiempo real. 


La eficiencia del equipo de centelleo líquido utilizado para 
la medición del tritio, disminuye debido al contenido de sales: 
- disueltas en las muestras. Para salvar dicho inconveniente se optó 
por destilar todas las muestras, lo que retrasó cerca de un mes los 
resultados, pero mejoró notoriamente la sensibilidad en 
las medidas. 


PARTE EXPERIMENTAL 
6.A — URANINA 
a) preparación de soluciones: 


Los estándares se prepararon con agua procedente de cada una 
de las transectas indicadas antes de la inyección de uranina. 


1) Solución madre: 


B 


Se disuelve 0.01009 de uranina llevando a volumen en matraz de 
1d 


2) Soluciones estandares: 
A partír de diluciones de la solución madre se preparan 
soluciones de concentraciones: 80,40 y 10 p/l (ppb). 


3) Solución blanco: 
Correspondiénte a agua de caúa una de las transectas extraída 
antes de la inyección de uranina. 


b) Calibración del equipo: 


Se efectuó siguiendo las indicaciones del manual para el caso 
en estudio. 

Se comprobó la linealidad en el rango de medida utilizando el 
ajuste de la curva de calibración que proporciona el software del 
instrumento. 


c) Análisis de las muestras 


Se fijaron las condiciones de lectura de modo que el valor 
registrado correspondiera al promedio de valores durante 30 seg a 
intervalos de 1 seg, realizado esto en forma automática por el 
instrumento. Esto evita errores potenciales de lecturas fluctuantes 
en el tiempo. 

Antes de iniciar las medidas se leen todas las soluciones 
estándares y blancos correspondientes, de modo úe verificar las 
condiciones. 

Luego se intercala la lectura de algún estándar cada 5 
lecturas de muestras; las cuales se leyeron por duplicado en 
diferentes días. 


DISCUSTON DE LOS RESULTADOS: 


1) Las muestras M1 y M2 no mostraron valores detectables por 
el equipo, sin embargo, la embarcación se encontro” encina de la 
mancha visible de fluorescencia. 

Según las características de la uranina, esta es perceptible a 
concentraciones de hasta 20 ppb, valor que supera el límite de 
cuantificación del equipo que es de 10 ppb. 

Esto muestra que al comienzo del experimento (al cabo de una 
hora)la uranina no está presente en la capa superior, por lo tanto 
el colorante tuvo que descender hasta capas más profundas, no 
habiendo tiempo suficiente para alcanzar la homogeneidad en la 
colunna de agua. Esto se comprueba con experiencias de laboratorio. 


2) Considerando los perfiles de conductividad de los días de 
muestreo se observa un salto marcado, sítuado en la mayoría de las 
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muestras a aproximadamente 3m de profundidad. Visto que el muestreo 
$e realiza a no más de 3,5m de profundidad (siendo la total entre 
5 y 10m), tomando una "columa" gue se llena con agua de esa 
profundidad, hace que las lecturas de uranina no Sean 
representativas del aporte del total de la distibución vertical del 
colorante. Si realmente se tomaran muestras en columna, los valores 
obtenidos se verian disminuidos. En efecto, el agua del estrato 
inferior es de elevada conductividad, lo que indica una entrada de 
agua oceánica libre de uranina. Tampoco es válido hacer un cálculo 
grosero de la cantidad de uranina en la franja analizada, ya que 
las muestras no son representativas del total del volumen de aguá 
en cuestión. Por lo tanto se recomienda hacer una variación en el 
proceso de toma de muestras si se desea realizarlas en columnas 


representativas. 


3) Tomando en cuenta las conclusiones anteriores se sugiere 
que la uranina viaja un una capa por encima del estrato de agua 
salada. Lo que está de acuerdo con un principio básico de la 
fisicoquímica que manifiesta que la solubilidaá de un compuesto 
orgánico disminuye apreciablemente el aumentar la concentración 
salina de las soluciones. 


4) Durante la fase de marea decreciente en el estuario, el 
flujo que circulaba hacia la zona donde se halla la Bahía de 
Montevideo encuentra una barrera, constituida por la escollera 
oeste (Ver fig.2). El flujo generalmente se divide en dos. Una 
parte se deflecta y continúa en dirección al este, paralelamente 
a la costa, mientras que la otra parte se deflecta penetrando hacia 
el interior de la Bahía, como muestra la fig. 

Esta parte del flujo arrastra los importantes aportes de efluentes 
industriales y urbanos que recibe diariamente la bahía y sale 
finalmente por la boca de acceso situada entre la escollera oeste 
y la escollera Sarandí. Durantre la fase de marea creciente 
generalmente se produce un desplazamiento global de masas de agua 
que probablemente sólo produce renovación en las másas de agua que 
se encuentran próximas a la boca de la bahía, sin al parecer el 


efecto descrito precedentemente durante la fase de márea 
descendente. 
En consecuencia la experiencia realizada con uranina en 


condiciones de bajamar, a pesar de los mencionados errores de toma 
de muestras, continúa siendo muy útil desde el punto de vista de 
la simulación del proceso de transporte de contaminantes aportados 
por la banía a la playas de Montevideo. 


5) Dado que la Uranina se destruye por la luz solar (longitides de 
onda corta) existiendo reportes de hasta un 50% de pérdida en 3 
horas de exposición (Feuerstein, D.L. Sellik, R.E. Fluorescent 
Tracers forn Dispersion Measurements, Journal of Sanitary 
Engineering, ASCE 89, 1-21 (1963) ), si ésta viajara en la 
superficie no hubiera sido posible su cuantificación varios días 
después de su inyección. Nuevamente se confirma la hipótesis de 
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que la uranina se desplaza en una capa profunda protegida de la luz 
solar por las capas superiores de agua. 


6) Naturalmente existen muchos materiales fluorescentes que 
tienen características similares a la uranina e interfieren con 
su medida. Pero cuando se utilizan filtros ópticos seleccionados 
cuidadosamente, se requerirían concentraciones más altas que las 
comunes y blancos fluorescentes más variables para alterar los 
resultados. En nuestro caso los blancos eran ceros al igual que 
muchas muestras, lo que descarta esta interferencia. 


7) En cuanto al PH sería recomendable asegurarse que sus 
valores naturales se encuentren entre 6 y 10, ya que a PH 5,5 cae 
abruptamente la lectura. 


8) Es de esperar resultados completamente diferentes con el 
tritio, debido a que al inyectarse como molécula de agua, facilita 
una rápida homogeneización en el sistema en forma de compuesto 
estable para los tiempos de experiencia utilizados. Con este 
compuesto se modelaría el comportamiento del agua de la Bahía de 
Montevideo, en cambio la uranina seguiría el curso de contaminantes 
de mayor peso molecular. 


6.B _ TRITIO 


La medición de la actividad de tritio en el agua tritiada 
(THO) se pensó efectuar en un comienzo sobre muestras filtradas en 
filtros de PVC (0,45pm) para eliminar el material en suspensión, 
a las que se le añadiría líquido de centelleo (picofluor Packard). 
Pero al investigador la relación entre eficiencia y conductividad 
(curva de calibración ) se halló una pronunciada inestabilidad en 
la zona de conductividades elevadas. A11í la eficiencia aumenta 
apreciablemente con el paso del tiempo (Ver fig). No obstante, se 
advierte una zona estable entre 0 y 10m8/cm. Para poder situarse en 
ese intervalo de conductividades se decidió destilar por columna 
cadá muestra. Se destiló la totalidad de las muestras obteniéndose 
conductividades comprendidas entre 15,5 S/cm y 8,76 mS/cm. Los 
viales preparados con picofluor fueron medidos y se comportaron en 
forma estable hasta siete días despúes de su preparación. Se 
intercalaron cada 25 viales, dos viales con estándares internos 
(patrones de THO) y dos viales con una solución "blanco" 
(destilada y con la concentración más baja encontrada en la 
medición de varias muestras tomadas antes de la inyección del 
trazador). 

Se obtuvieron los resultados que muestra la tabla 1I adjunta 
(ver Apendice 3). 


1—RADTOPROTECCION 


dl. 


En este tipo de experiencias, donde se utilizan trazadores 
radiactivos para el estudio y evaluación de problemas 
hidrodinánicos y ambientales, es imperioso justifacar plenamente la 
metodología empleada. A tales efectos se elaboro un informe en 
conjunto con la Divigión de Radioprotección de la Dirección 
Nacional de Tecnología Nuclear, donde se evaluaron las ventajas y 
desventajas del método. En este caso la evaluación demostró que 
la tarea se encuentra plenamente justificada. (Ver Apendice 1) 


8-MODELIZACION Y SIMULACION DIGITAL DEL SISTEMA 


La planificación de los experimentos con  trazadores 
radioactivos involucró una fase prevía de modelización y simulación 
digital de los procesos de transporte involucrados. Se tuvo en 
cuenta la estratificación vertical, la acción de los vientos, las 
profundidades y otras condiciones de la zona, todo ésto en el 
contexto de un nodelo en el que se define un piso difusivo para el 
transporte del trazador y luego se integra en profundidad. Se 
emplea un tensor de coeficientes de dispersión en el que aparecen 
parámetros relacionados con procesos advectivos no resueltos en la 
escala del modelo del sistema, tales como corrientes hacia y desde 
la costa (variables en profundidad y asociadas a los vientos). 
Los efectos debidos a las modificaciones en la profundidad se 
integraron en un término advectivo adicional. [(Suárez, Badano, 
Dellepere y otros,1993)]. 

Para estimar la posición del piso difusivo se enplearon 
correlaciones clásicas entre el coeficiente de dispersividad 
vertical y el número de Richardson, calculando valores de 
orientación para este último en función de la profundidad, en base 
a estimaciones locales de la frecuencia de Brunt-Vaisala a partir 
de perfiles de salinidad y temperatura, y de valores de operación 
de la frecuencia de Kelvin-Helmholtz obtenida a partir de 
experiencias con biplanos. 


9-CONCLUSTONES Y RECOMENDACIONES 


El agua de la Bahía de Montevideo circula a lo largo de una 
franja costera cuya anchura puede oscilar entre 1 y más de 5 kms. 
, dependiendo de la zona, con la zona más estrecha situada a la 
altura de la saliente de Punta Carretas. 

Los procesos de transporte son muy rápidos a lo largo de la 
costa: en dos días una partícula que parte próxima a la boca de la 
Bahía puede alcanzar la Playa Carrasco, a más de 15kms.. úe 
distancia. 

La primera detección de tritio se efectuó a las 2 horas de 
ocurrida la pleamar siguiente a la inyección del trazador. En ese 
entonces, el agua marcada había llegado ya a la Playa Ramírez 
(transecta 1) con una concentración alta en un punto próximo a la 
costa. Al día siguiente, no se encuentra tritio en la transecta 2 
(Playa Pocitos) ni en Ja zona más próxima a la costa de la 
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transecta 3. Pero se halla el trazador en la zona más alejada de la 
trangecta 3 y en la transectal. 

Como consecuencia, el efecto de las descargas locales úde 
contaminantes es relativamente más importante allí respecto a otros 
puntos de la costa. 

Los resultados de las medidas de conductividad del agua 
corroboran esta conclusión. En efecto, durante los días 24 al 28 de 
marzo de 1993, se produjo una invasión de agua de baja salinidad 
proveniente del Oeste. Durante ese tiempo, la conductividad en la 
Playa Pocitos no descendió a los valores tan bajos medidos en la 
Playa Ramírez (4 mS/cm) situada, respecto del flujo proveniente de 
Oeste, antes de la saliente de Punta Carretas. 

Á partir del día 28 se produjo una invasión de agua más 
salina, con conductividades mayores que 35 mS/cn. 

La aparición de masas de agua de conductividades diferentes 
marchando en sentido opuesto y mezclándose con el agua marcada no 
llegan a producir una renovación apreciable del agua de la zona 
costera investigada, puesto que al final de esos procesos, el 
tritio siguió estando distribuído a lo largo de la costa, con 
concentraciones variables, más pequeñas, pero medibles sin 
dificultad. 

Al parecer entonces, los fenómenos migratorios de la interfase 
entre agua dulce y agua salada que se producenen la zona costera 
del Departamento de Montevideo se caracterizan por un mezclado 
relativamente lento, con tiempos de permanencia comprendidos entre 
5 y 10 días. 

En suma, los resultados de la marcación con tritio del agua 
que sale de la Bahía indican claramente que, en las condiciones 
existentes durante la campaña, la mayor parte de la contaminación 
gue se observa en las playas se debe al aporte del agua contaminada 
de la bahía. Además hay un cierto aporte debido a los vertidos que 
se efectúan en las mismas, en las zonas de agua relativamente 
atrapada. 

La influencia que sobre la contaminación podría ejercer el 
emisario subacuático de Punta Carretas, en las playas de 
Montevideo, al parecer es muy pequeña, incuso despreciable en las 
condiciones en las que se efectuaron las experiencias. No obstante, 
son necesarios nuevos estudios para definir este punto y extender 
la apiicabilidaá de la conclusión referente al rol de la Bahía en 
la contaminación de las playas. 


Como consecuencia de este estudio se desprende que la escala de 
tiempo asociaáa a los procesos advectivo-dispersivofla lo largo de 
la costa del Departamento de Montevideo, desde la boca de la bahía 
hasta la playa Carrasco, es del orden de tres días. 

Un elemento, que eventualmente sufra degradación (o muerte, en caso 
de microorganismos) pero cuya escala de tiempo sea de orden 
superior a la del transporte lo consideraremos como un elemento 
consérvativo, a los efectos de esta línea de investigación. Como 
consecuencia un elemento de este tipo puede alcanzar el límite de 
las aguas costeras del Departamento de Montevideo, con un 
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decaimiento poco significativo. 

En vista de la importancia del las conclusiones de este estudio, se 
recomienda continuar los trabajos enprendidos en dos direcciones. 
Desde el punto de vista teórico tendría interés el simular la 
evolución de la concentración promedio de un contaminante en la 
bahía, calculando la cantidad que es expulsada durante cada ciclo 
ás marea. Esto podría emplearse como entrada a un modelo que 
representa el proceso de advección-dispersión a lo largo de la 
costa y permita estimar concentraciones para otros escenarios 
diferentes al que sirvió de marco a la campaña experimental que 
describimos. 

Para ello es necesario prestar particular atención a los fenómenos 
de transporte en la boca de la bahía . Se puede obtener información 
valiosa , adicional, previa al desarrollo del modelo del transporte 
hacia y desde la bahía, a partir de mediciones con trazadores, 
biplanos y  correntógrafos. Para ajustar parámetros pueden 
emplearse, además, los resultados experimentales obtenidos en la 
campaña de 1991 (Plata y otros,1992) 

Desde el punto de vista experimental convendría estudiar el 
comportamiento de las aguas residuales vertidas en la bahía por los 
arroyos Miguelete y Pantanoso . 

Sería conveniente marcar con trazadores la desembocadura de estos 
arroyos. Log resultados obtenidos podrían compararse con los 
resultados que arroje una eventual simulación digital, basada en 
un software adecuado. La atención podría dirigirse entonces hacia 
la determinación de la función de influencia que relaciona un 
vértido de una sustancia en la desembocadura de los mencionados 
arroyos con la posterior presencia de dicha sustancia en la 
desembocadura de la bahía 
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INFORME PRELIMINAR DE SEGURIDAD RADIOLOGÍCA PARA LA CAMPAÑA DE 
HIDROLOGÍA 1993 
DNEN 


1. INTRODUCCION 
IT. OBJETIVOS 


111. JUSTIFICACIOÓN 
a. Procedimiento con tritio 
b. Procedimiento con 1-131 
c. Ventajas del procedimiento empleado 


1V. METODOS DE CAMPO 
a. Marcado con tritio 
b. Marcado con 1-131 
c. Medidas de seguridad y monitoreo ambiental 
ád. Medidas de seguridad en situación no planificada 


V. CONCLUSION 


T. INTRODUCCION 


La Bahía de la ciudad de Montevideo, como muestran los datos 
de la campaña anterior realizada por la DNTN (del 2 al 13 de 
diciembre de 1991), cuenta con un volumen medio total de agua de 
48.61 millones de metros cúbicos para una superficie de 15.92 
millones de metros cuadrados, lo que da una profundidad media de 
3,05 metros. 


En «lla desembocan los arroyos Pantanoso y Miguelete, y 
descargan los colectores: 


- el interceptor de 3370 m de extensión, que bordea por el 
norte y el oeste de la ciudad vieja. Nace en las calles Juncal y 25 
de Agosto y termina en la Rambla Sur y la calle Guaraní. 

- él colector construido por la Avenida Rondeau y la calle 
Zelmar Michelini, de 3300 m de extensión (de la cual 1416 m son en 
túnel). Nace en la zona de Arroyo Seco y descargá en la Rambla Sur 
y la calle Paraguay, 

- y un vertimiento en la zona oeste de la industria Ancap (ver 
figura 1). 


Referente a la contaminación que provocan estas descargas 
puntuales, cabe destacar lo siguiente: 


1) Arroyos Pantanoso y Miguelete: dadas las características 
morfológicas y físicas de las cuencas de ambos arroyos, 
el régimen hidrológico es insuficiente para soportar las 
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cargas aportadas a dichos cursos. Las fuentes principales 
de contaminación son las industrias, las descargas 
directas Óó indirectas de saneamiento y el aporte de 
residuos sólidos producidos y manejados en la cuenca. 

2) si bien ambos colectores tienen sus puntos de vertimiento 
fuera de la bahía, existen a lo largo de los colectores 
una serie de aliviaderos con la finalidad de evacuar los 
caudales que superen los caudales de diseño. Debido al 
gran aumento de los caudales que llegan al colector en 
los últimos años (mayor población servida) el 70% del 
agua servida esta vertiendo continuamente en la bahía 
(por los aliviaderos) y apenas llega á la costa sur el 
30% del caudal de tiempo seco. 

3) Enisaríio Ancap: debido a lo antes mencionado, y sumado a 
la actividad portuaria existente en la bahía, explica la 
alta contaminación orgánica, patógena y tóxica actual de 
la misma. 


Cabe señalar que considerada en su conjunto, el agua de la 
bahía se renueva, en condiciones meteorológicas normales, con un 
tiempo de renovación nominal (relación volumen/caudal) que debe 
oscilar, según las circunstancias, entre 15 y 40 horas (según datos 
de la anterior campaña). 


TI. OBJETIVOS 
Los objetivos planteados para esta segunda campaña son: 


1) A pedido de la Facultad de Ingeniería, el estudio de la 
posible influencia de la contaminación de las aguas de la 
bahía sobre las playas de Montevideo. 


2) A pedido de la Intendencia de Montevideo el estudio de la 
pluma que emite el emisario de Punta Carretas. 


Es importante destacar que ambos estudios son de gran 
relevancia nacional y que hasta el presente los estudios realizados 
por otros organismos no han sido consistentes en sus apreciaciones 
finales. 


111. JUSTIFICACIÓN 


En función de los dos objetivos planteados y teniendo en 
cuenta los sistemas que deben ser estudiados, la metodología con 
tragadores radiactivos es la mas aconsejable para obtener la 
información de naturaleza global requerida. Los criterios 
necesarios para la justificación de la práctica se detallan a 
continuación. 


318 


á. Procedimiento con tritio 


Como se ha expuesto anteriormente la bahía de Montevideo 
recibe diariamente importantes caudales de liquidos cloacales y de 
uso industrial, lo que explica una alta contaminación en ella. En 
momentos de marea baja cuando ingresa agua a la bahía, comienza a 
aumentar su volumen acumulándose en ella los vertimientos 
mencionados. En momentos de marea alta esta masa de agua 
contaminada sale de la bahía. 


Hasta el momento los estudios existentes no permiten obtener 
información sobre que tipo de influencia genera el agua de la bahía 
sobre las playas de Montevideo. 


La metodología propuesta para realizar el estudio pedido, 
implica el vertimiento de 100 Ci de agua tritiada en la boca de la 
banía, obtienéndose al cabo de minutos concentraciones aámisibles 
(ver calculo en Métodos de campo). 


Cabe mencionar como ejemplo el caso del arroyo Pantanoso: el 
siguiente cuadro compara los valores de algunos parámetros de 
contaminación, muestreados por la DI.NA.M.A., en puntos de muestreo 
cercanos a la báhía, con respecto a los límites de los parámetros 
de contaminación correspondiente a la clase 3 del Código de Aguas. 


EA 
9.025 
31.73 3.7805 


Situación semejante ocurre con el Arroyo Miguelete, cuyo 
régimen hidrológico es insuficiente para soportar la carga aportada 
por aproximadamente 494000 habitantes equivalentes. Solo a los 
efectos de comparación, usando el valor guía manejado en estudios 
recomendado para sostener el uso como cuerpo receptor de residuos 
crudos estaría dando que la cuenca del Miguélete tendría solo 
capacidad suficiente para recibir la carga de 15000 habitantes 
equivalentes. 


PARAMETRO 


PTO. MUESTREO 
Piuvial Alaska 


) 


DBO(mg/L) | CF(/100mg) 


55000 


856 


LÍMITE MAX. 


AUMENTO 


con referencia a la situación del Arroyo Pantanoso se han 
alegido solo 4 parámetros de contaminación a los simples efectos de 
mostrar el gran aumento entre Jas medidas tomadas y los máximos 


permitidos; pero esta situación se repite para los demás 
parámetros. 


Se podría pensar en la utilización de otro tipo de trazador 
para realizar el estudio, por ejemplo los trazadores químicos como 


A 


ser dicromato potásico, uranina o rodamina, sin embargo presentan 
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algunos inconvenientes que hacen muy dificultosa su utilización: 
Para cualquiera de estos trazadores se necesitarían cantidades que 
oscilarían en unos 500 Kg (A. Plata Bedmar, "Isótopos en 
hidrología", Ed. Alhambra, 1972), tomando el precio mas bajo que 
corresponde a la uranina = US$ 77/Kg el estudio tendría un costo 
del oráen de US$ 40000, solo por la adquisición de estos productos. 
El costo de 100 Ci de tritio es de US$ 9000; además de contar con 
una extraordinaria sensibilidad de medida. 


b. Procedimiento con 11-131 


En este caso, se utilizará un radioisótopo con bajo valor del 
período de semidesintegración del 1-131 (8 días), con una actividad 
total a inyectar de 6 Ci. Cabe destacar la ventaja que posee este 
trazador frente a otros, al ser un emisor gamma, lo que permite 
topografiar la pluma in situ. 


c. Ventajas del procedimiento empleado 


1) En el caso del 11-131, a) posibilidad de medida in situ, b) 
viáa limitada del trazador y c) posibilidad de emplearse en 
diferentes formas químicas sin que se modifique la 
sensibilidad de la medida. Esto es importante, con vistas a 
conseguir que el comportamiento del trazádor se asemeje al de 
ias moléculas de agua. 


2) En el caso del tritio, al ser'agua tritiada tendrá el mismo 
comportamiento que el agua a estudiar. 


IV. METODOS DÉ CAMPO 


En ambas operaciones de campo, un oficial de radioprotección 
supervisará la tarea. con monitores apropiados in situ. Se 
utilizarán detectores personales para medición de dosis por 
irradiación externa. Se contará con equipos de seguridad (guantes, 
soluciones decontamidoras). 


a. Marcádo con tritio 


Se procederá al marcado con tritio mediante una inyección 
puntual localizada a 200 m al sur de la escollera Sarandí y a 200 
m del extremo oeste desde una barca. 5e procederá al rompimiento de 
la ampolla conteniendo el agua tritiada 2.5 m debajo de la 
superficie del agua con un dispositivo especialmente construído 
para ese fin (probado en frío con anterioridad). La barca 
permanecerá como máximo 10 minutos en el lugar de vertido. Antes 
del vertido de la actividad, se comprobará mediante vertido de 
uúuranina y mediante flotadores, el sentido del flujo de las 
corrientes. Teniendo en cuenta que el límite anual para ingestión 
de tritio para seres humanos (ALI) es de 32.4 mCi, de acuerdo a las 
recomendaciones del Technical Report n*%324 del OIEA y de las nuevas 


20 


recomendaciones del ICRP n*60, y considerando una dilución en un 
cubo de 3 m dé lado, se obtiene una concentración de 3.7 mCi/L. 
Esta hipótesis es altamente conservadora, del lado de la seguridad. 
Para llegar a los límites recomendados, habría que ingerir 3 litros 
de esa agua. Cabe resaltar que el agua tritiada tiene un 
comportamiento químico idéntico al del agua no tritiada, para las 
condiciones de interés. 


Simultáneamente, se hará un trazado con uranina lo cuál 
permitira preáécir y visualizar la futura evolución de la nube de 
tritio. Se verificará que en la zona de la inyección no halia otras 
embarcaciones. 


b. Marcado con 1-131 


Se realizarán 3 inyecciones de 2 Ci cada una al cabo de un 
día. Cabe señalar que los difusores se encuentra a 4 m de 
profundidad y a 2.9 km de la costa. 


c. Kedidas de seguridad y monitoreo ambiental 


En cuanto al transporte y almacenamiento de las fuentes, las 
sustancias radiactivas estarán en todo momento bajo el control de 
un oficial de radioprotección asignado, única persona capaz de 
autorizar y supervirzar todos sus movimientos. 


En relación al impacto de la práctica en las playas de 
Montevideo, sé efectuará un muestreo del agua de todas las playas 
de Montevideo para evaluar la actividad especifica de tritio e iodo 
en ellas, antes y después de los vertimientos de material 
radiactivo. Los análisis de tritio se efectuarán en la cátedra de 
Radioquiímica, Facultad de Química, coordinados por el Laboratorio 
de Radloanálisis, y los de jodo en este último. 


á. Medidas de seguridad en situación no planificada 


La comprobación por parte de los operadores de campo de una 
situación no planificada que involucre riesgo potencial de 
exposición o contaminación, deberá ser comunicada de inmediato al 
responsable del proyecto (Prof. R. Suarez) y al oficial de 
radioprotección asignado por medio de sistemas de comunicación 
adecuados, que deberán chequearse al comienzo de la práctica. Ambas 
personas tomarán de mutuo acuerdo las decisiones del caso, para lo 
cual se presenta una lista de casos a considerar: 


o contaminación con agua tritiada: el operador deberá de 
inmediato descontaminarse con agua no tritiada, e ingerir 
abundante cantidad de agua no tritiada, suplementados con 
sustancias de efecto diurético. Cuanto antes se realicen estos 
procedimientos, mucho menor será la dosis producida. 

- contaminación con diste deberá procederse a la 


descontaminación con abundante agua jJabonosa, y Sl se 
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sustificara, proceder al bloqueo de la glándula tiroides con 
1 no radiactivo. 


V. CONCLUSION 


El presente informe ha sido confeccionado en base a la 
interacción entre los responsables del proyecto y la División 
Protección Radiológica y Seguridad Nuciear. En él se busca 
establecer una base para la aprobación desde el punto de vista de 
la seguridad radiológica de las prácticas mencionadas. La idea 
inicial de la auditoría ha sido contemplada ya que en todo momento 
el interés fue el intercambio de ideas y el uso del sentido común 
para mejorar la seguridad. Se contemplaron básicamente tres 
aspectos claves: seguridad de las sustancias radiactivas, seguridad 
de los operadores, y seguridad del público y medio ambiente. Lo 
anteriormente expuesto muestra que, en las condiciones de seguridad 
radiológicas planificadas (descritas en este informe), la práctica 
puede ser efectuada de acuerdo a los criterios de seguridad 
aceptados. 


Han participado en la preparación de este informe, 


por la División Protección Radiológica y Soquridad Nuclear de DNTN: 
»>A4ido Badano 
por el grupo de Hidrología de la DNTN: .«=pRoberto Suarez 
->»Alejandro Dellepere 
>Martin Barreiro 
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Tabla 1, Resultados de las concentraciones de trítio obtenidas 
para las muestras tomadas después de la inyección de este 
trazador. Análisis efectuados en la Facultad de Ciencías 
Quimicas de la Universidad de Montevideo. Las concentraciones 
vlenen dadas en Bq/l. 


PUNTO 23/3/93 24/3/9324/3/93 25/3/93 26/3/93 28/3/93 1/4/93 
mafíana tarde 


10 44,5 

11 9,48 4,99 7,27 19,46 13,30 3,57 25,37 
12 0,00 4,37 4,37 13,67 3,69 7,88 30,54 
13 4,93 9,30 11,45 13,11 6,40 8,80 35,22 
14 3,26 19,15 7,45 12,36 4,00 7,76 31,65 
15 7,20 17,24 6,10 5,05 10,22 25,06 
16 4,93 0,00 5,05 

21 7,33 5,42 4,86 10,59 8,80 6,26 27,21 
22 0,00 0,00 7,76 6,31 9,85 9,85 27,77 
33 0,00 0,00 11,21 17,55 8,00 10,47 27,15 
24 5,54 0,00 0,00 12,13 6,46 10,04 29,74 
25 0,00 10,71 3,88 8,80 24,81 
26 10,71 3,69 

31 0,00 14,47 6,53 21,96 
32 4,31 8,80 6,59 4,37 25,43 
33 5,91 16,62 6,90 7,27 23,69 
34 10,53 18,10 7,82 3,02 23,63 
35 12,15 11,64 5,42 3,69 29,55 
36 17,12 19,09 

40 7,45 

41 6,90 14,65 10,65 27,46 
42 4,49 10,59 6,03 8,56 9,17 
43 0,00 11,70 7.63 6,90 31,28 
54 4,56 14,90 9,05 0,00 23,46 
45 10,47 8,74 6,96 0,00 18,22 
46 5,23 8,62 
51 11,51 58,00 5,29 38,73 
52 11,70 12,19 9,30 23,46 
53 9,24 9,87 4,06 

54 7,88 5,36 13,73 21,86 
55 9,67 6,34 4,62 34,29 
56 15,21 8,43 5,23 26,17 
57 8,58 1,37 33,83 
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ESTUDIOS SEDIMENTOLOGICOS EN LA REGIÓN DE LA BAHÍA 
DE MONTEVIDEO, CON EL EMPLEO DE TECNICAS NUCLEARES 
Y TRAZADORES RADIACTIVOS Y FLUORESCENTES 


Geraldo Godinho Pinto”, Jefferson Vianna Bandeira”, Roberto $. Antola?, 
Alejandro Malek”, Beatriz S. Amei enda”, Maria Luisa B. Parada? 
Hugo A. Fraga?, Diego Del Monte! Y Anna Firpo”, Walter Cabral? 

R. Odino?, Mary Postiglione*”, Graciela De León Magnani”, Aldo Cruz”, 

Jaime B. Gorfain”, Heber Longoni”, Hector Cane', J orge D'Angelo”, 


Felipe Ramirez”, Luis A. Tapia'”, Silvia M. Tasende”, Pablo Garcial* 


INTRODUCCION 

Este trabajo está relacionado con los estudios sedimentológicos con la 
utilización de trazadores radiactivos: 'Au (oro-198) y Te (tecnecio 99) y 
trazador fluorescente (fluoresceína sódica) en la Bahía de Montevideo, 
incluyendo la región del Río de la Plata cercana a ella (Figura la), hecho en el 
marco del proyecto de cooperación técnica del Organismo Internacional de 
Energía Atómica (OIEA) URU/8/009, con los objetivos de evaluar el 
comportamiento del sedimento fino dragado en el puerto de Montevideo y 
vaciado aguas afuera de la bahía, en el Río de la Plata, y la dinámica del agua y 
del material en suspensión en el interior de la misma, siendo la dinámica del 
agua estudiada utilizando trazador fluorescente. 

Estos estudios también contemplaron el empleo de sonda para medición 
de la densidad de limo del fondo, basada en la transmisión de la radiación y. La 
sonda fue desarrollada en el CDTN/CNEN-Brasil, por el autor principal del 
presente trabajo. Como soporte a estos estudios fueron realizadas mediciones 
de corriente y de parámetros del agua: concentración in situ de sólidos en 
suspensión, conductividad, temperatura y salinidad. 


A A a 

(1) - CDTN/CNEN - CP 941 - 30123-970 Belo Horizonte, MG, Brasil. Tel: 031 499-3127; 
Fax: 031 499-3249. E-mail: ¡vbCurano.cdtn.br 

(2) - Dirección Nacional de Tecnología Nuclear (DINATEN) — Mercedes 1041, CP 11.100 
Montevideo, Uruguay. Tel: (+598 2) 9026919; Fax: (+598 2) 9021619 

(3) - Ministerio de Transporte y Obras Publicas (MTOP) — Dirección Nacional de Hidrografía 
Rincón 575, 2” Piso, Montevideo, Uruguay. Tel: (+598 2) 9064783; 
Fax: (+598 2) 9064667 

(4) - Administración Nacional de Puertos (ANP) - División de Operaciones Marítimas 

(5) - Centro de Investigaciones Nucleares (CIN) — Universidad de la República 

(6) - Universidad Católica del Uruguay 
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DESCRIPCION DE LA REGION 

La Bahía de Montevideo es pequeña (16 km?) y esta situada en la orilla 
izquierda y en el tramo inferior del Río de la Plata (Figura 1b) adonde hay una 
gran variación de salinidad (Sihn, 1989; Plata et alli., 1992). La circulación del 
agua en la bahía es sobretodo en el sentido horario, siendo la profundidad 
natural en el interior de la bahía de 1 a 2 m, teniendo el limo como material 
predominante. 

El acceso al puerto es hecho por un canal dragado a 10 m de profundidad 
y con longitud de 10 km. En el interior de la bahía este canal sigue hacia 
Noroeste para permitir el acceso de buques a las instalaciones petroleras de 
ANCAP (Figura la). 

El sitio de vaciado de dragado es en el Río de la Plata, al Este del canal 
de navegación, cerca de 4 km desde la extremidad sur de la escollera Oeste 
(Figura la). El material dragado es sobretodo limo y arcilla, siendo que el 
fondo presenta una capa extensa con variación de densidad para el estudio de la 
cual fue desarrollada la sonda nuclear. 

Debido a la existencia de industrias en la región de la bahía, hay un 
vertido de materiales contaminantes directamente en la misma y sobretodo a 
través de los ríos Pantanoso y Miguelete. Los contaminantes son transportados 
por el agua y también en los sedimentos finos en suspensión (sobretodo 
metales pesados), debido a la adsorción. Así, el estudio de la dinámica del agua 
y de los sedimentos en suspensión son importantes para los estudios 
ambientales en la bahía. El proyecto del OJEA propició la transferencia de 
tecnología en el empleo de técnicas nucleares para los estudios mencionados. 


ESTUDIOS REALIZADOS 

Para el estudio del comportamiento del vaciado de dragado en el Río de 
la Plata fueron hechos experimentos con el empleo de 198 Au marcando la carga 
de la Draga 7 de la Administración Nacional de Puertos (ANP), que es Una 
draga del tipo TSHD (“Trailing Suction Hopper Dredger”) con 2.690 m” de 
capacidad en la tolva. La marcación consistió en promover la adsorción del 
'*Au en un volumen de 50 litros de limo homogeneizado en uno circuito 
cerrado accionado por una bomba hidráulica, y mezclado con el trazador 
empleando una segunda bomba flujo controlado, justo antes de ser vertido en la 
tubería de alimentación de la tolva de la draga durante su cargamento (Tola, 
1982). Fueron hechos dos marcados de limo en la Draga 7 (Bandeira, Pinto, 
1996), siendo el material vertido en el punto de inyección PI (Figura la). 

El estudio del comportamiento en suspensión del material vaciado fue 
hecho con el empleo de contadores de centelleo a 4 y 6 m de profundidad, 


315 


remolcados por un barco ubicado electrónicamente y, en el fondo, después de 
la detección en suspensión, con un contador de centelleo acoplado a un trineo. 

Para el estudio del comportamiento del limo en suspensión existente en la 
bahía, fue empleado el mismo procedimiento de marcado para la draga pero 
haciendo una inyección sub-superficial del sedimento marcado, con una 
concentración cerca de la natural de los sedimentos en suspensión en la bahía, 
para no crear corrientes de densidad. La evolución de la nube de material 
marcado fue rastreada con el empleo de un contador de centelleo, ahora a 0,5 m 
bajo de la Superficie (Bandeira, Pinto, 1997). Debido a la imposibilidad de 
recibir el 'Au en tiempo hábil para realizar los trabajos, fue empleado, en 
experimento inédito, el PTc en el marcado de sedimento fino. La idea del 
empleo de Este trazador fue de la Química Beatriz Souto. Los estudios de 
fijación del Tc al sedimento fueron hechos en los laboratorios de DINATEN. 

La fluoresceína sódica fue mezclada al limo marcado con “Tc para se 
tener una inyección simultánea de los dos trazadores. Para la detección del 
trazador fluorescente, cuya evolución representa la circulación del agua en la 
bahía, fue empleado un fluorimetro Turner, modelo 10-AU, teniendo la toma 
del agua también ubicada a 0,5 m de la superficie. Flotantes fueron empleados 
para facilitar la visualización del desplazamiento de la nube (Bandeira, Pinto, 
1997). 

Anteriormente a los trabajos arriba nciicionados, fueron hechas 
mediciones de corriente (perfiles verticales), conductividad, temperatura, 
salinidad y concentración in situ de sedimentos en suspensión, utilizando para 
este parámetro, un turbidímetro óptico. 

Para la determinación del perfil de densidad de limo depositado en el 
fondo de la bahía, fue desarrollada una sonda basada en el principio de la 
transmisión de la radiación y (Figuras 2a y 2b), presentando las siguientes 
características: 

e. Material: latón y acero inoxidable; 

e Fuente: Am (americio) Amersham tipo X101, código AMC.P4 (14 mCi) 
con posibilidad de utilización de fuentes de 10 a 100 mCi; 

e Distancia fuente-detector: variable de 8 a 30 cm; 

e Haz de radiación con posibilidad de ajuste en los tres ejes; 

e Detector de radiación: Minekin, modelo 9302; 

e Sonda de presión Eclipse (0 a 10 bars); 

e Rango de mediciones: 1,0 a 1,3 ton/m'; 

+ Peso del equipo: 11,5 Kg. 


RESULTADOS 
Experimentos con A y en la Draga 7 
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En el primero experimento se obtuvieron: 

1) el limo marcado en suspensión fue desplazado para SO, de acuerdo con la 
dirección de las corrientes más profundas medidas antes de la inyección; 

2) debido al nivel de apertura de las puertas de la draga (4 a 5 m), respecto a la 
profundidad total del sitio de vaciado (8,3 m), y también bajo la influencia del 
momento vertical descendente y a la corriente de densidad producida por el 
vaciado, una gran parte del sedimento se desplazó rápidamente para la región 
más profunda y se depositó en el fondo cerca del sitio de vaciado, lo que fue 
verificado por una detección de fondo, justo después de la detección en 
suspensión. Esto no significa que el material de vaciado permanezca en el 
fondo. 

El segundo experimento no fue concluido debido a condiciones 
climáticas afuera de la bahía (fuerte corriente de crecida asociada a fuertes 
vientos del cuadrante sur, con intensidades alcanzando 70 a 75 km/h). Estas 
condiciones impedían las maniobras necesarias a una buena navegación de 
detección. Así, una vez que el limo en la tolva de la Draga 7 ya fuera marcado 
en el ante puerto, fue decidido vaciarlo y volver al puerto sin hacer la 
detección. 

Experimentos con Te 

En la búsqueda del substituto para el 'Au, fueron desarrollados en el 
poco tiempo disponible en el cronograma, estudios preliminares de marcación 
del limo con “Tc. Los rendimientos de fijación obtenidos en dos experimentos 
utilizando dos reductores disponibles fueron: 


(1) Na2S>03 mezclado con NH¿OH.: 54,0%; (pH > 9) 

(2) SnCl> : 71,4%. 

Entonces, para el marcado del limo, la reducción del Tc con SnCl, fue 
escogida. 

Dadas las condiciones excepcionales del presente trabajo y a la media 
vida corta del “Te (6,0 horas), la actividad inyectada fue solamente de 45 mCi. 
Así mismo, fue posible acompañar la evolución del trazador por cerca de 
treinta cruces transversales de la nube, lo que no habría sido posible si el 
estudio hubiera sido de vaciado de dragado pues, bajo la influencia del 
momento vertical descendente y a la corriente de densidad producida por el 
vaciado de una draga “hopper” o un ganguile, la mayoría del sedimento se 
desplazaría rápidamente para la región más profunda, como sucedió en el 
vaciado hecho con la Draga—7 de la ANP. 

Los conteos en los picos mayores, durante la detección, fueron del orden 
de 20 veces el blanco: de 450 a 500 cuentas por 10 segundos (cp10s), mientras 
el blanco era del orden de 15 a 25 cp10s. Si la actividad inicial de “Te tuviera 
sido de 1200 mCi, actividad normalmente disponible para aplicaciones médicas 
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todos los lunes a las 1U:UUH en Montevideo, y la inyección hecha en este 
mismo momento, los conteos en los picos podrían tener alcanzado el orden de 
las 13.000 cp10s, lo que representaría cerca de 650 veces el blanco. Así se abre 
una gran posibilidad de empleo del Tc en estudios sedimentológicos en la 
Bahía de Montevideo y otros sitios, con la condición de que se tenga una buena 
fijación del trazador en el limo. Por lo tanto, estudios de fijación y del 
comportamiento físico-químico y hidrodinámico del limo marcado deben de 
proseguir en laboratorio. 

Detección de fluoresceína sódica 

Fueron hechos cerca de cincuenta cruces de la nube y obtenidos cerca de 
treinta picos bien definidos. El blanco medido fue del orden de 656 + 18 y 
fueron obtenidos picks hasta 1000. Los picks, corregidos del tiempo de 
bombeo, fueron simultáneos con aquellos debido al trazador radiactivo. Fueron 
tomadas también posiciones sucesivas de los flotantes a lo largo de la 
detección. La nube tuvo un dirección general hacia Norte. 

En conclusión, los resultados de los trabajos con el empleo de limo 
marcado con Te y el agua con trazador fluorescente permitieron seguir la 
circulación del agua y de los sedimentos en suspensión desde la apertura Oeste 
de la bahía hasta la región vecina al sur de la Isla Libertad (Figura 1 a). 
Empleo de la sonda nuclear 

El equipo fue calibrado en los laboratorios de DINATEN en el rango de 
densidades (1,0229 a 1,1869 ton/m?), con sedimento tomado del fondo de la 
bahía, en la región del puerto. Los valores R¡/Ry adonde R¡ y Ro en (cpm) en el 
medio y en el agua, respectivamente, se relacionan con la densidad del medio a 
través de la ecuación: : 


R¡/ Ro= expGpUpx) (1) 
adonde: 
u” = coeficiente de absorción másico del medio (cm”/g) 
p = masa específica de la mezcla (g/cm' ) 
x = distancia fuente-detector (cm) (x empleado = 11,5 cm) 
La ecuación de regresión lineal para los valores experimentales es: 
In(R¡/Ro) = 1,9002 - 1,8679p (2) 


con el coeficiente de correlación r = -0,9702. 

Las figuras 3a y 3b presentan los perfiles de densidad obtenidos en la 
Dársena II del puerto (Figura la), apuntando un espesor de cerca de 4,0 m entre 
las densidades de 1,02 y 1,20 ton/m? de la mezcla de limo y agua, 
respectivamente, las densidades indicadas por el eco sonda de 200 KHz 
(profundidad náutica) y la profundidad de maniobra sin riesgo para los buques. 
La utilización de la sonda en distintos puntos de la región portuaria permitirá 
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dibujar la curva de densidad = 1,20 ton/m', lo que implica en un gran ahorro en 
términos de gastos en dragado. 
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Abstract - This work is related to the sedimentological studies with the 
application of radioactive tracers, carried out in the region of Puerto Carmelo, 
Río de La Plata - Uruguay, under the TC Project of the International Atomic 
Energy Agency (IAEA) URU/8/009, with the objective of evaluating the bed 
load transport of sand, in the region where a new access channel to the harbour 
of Carmelo, was designed. These studies, performed together hydraulic 
measurements, were important for the comprehension of the behaviour of the 
sediment transportation under the action of the hydrodynamic agents (currents 
and waves); the determination of the bed load rate and ¡ts application to the: 
evaluation of the maintenance dredging; the definition of the optimum path for 
the access channel, based on the deposition process of the injected tracer, 30 
cm above the bottom. The maximum accretion rate for the designed channel is 
around 11.000 m'/year. 


Palabras llave (key words) - trazador radiactivo (radioactive tracer); 
sedimento (sediment); dragado (dredging); azolvamiento (accretion); Río de la 
Plata. 
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